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EXECUTIVE  SUMMARY 


We  are  evaluating  the -extent  to  which  biodiversity  in  harvested  areas  is  affected  by  the  size  and  number 
of  tree  clumps  that  are  left  standing  after  harvest.  Twelve  experimental  cut-blocks  were  created  in  each 
of  two  study  areas.  One  experimental  area  is  in  pine-dominated  forests  near  Whitecourt,  Alberta,  and  the 
other  is  in  aspen-dominated  forests  near  Grand  Prairie,  Alberta.  Within  each  of  the  cut-blocks  in  the 
experimental  areas,  3%  of  the  pre-harvest  merchantable  trees  and  snags  were  retained  in  clumps.  One 
third  of  the  experimental  cut-blocks  had  clumps  of  0.60-0.75  ha  (1.48-1.85  acres),  one  third  had  clumps 
of  0.14-0.15  ha  (0.35-0.37  acres),  and  one  third  had  clumps  of  0.03-0.04  ha  (0.07-0.10  acres).  This 
research  complements  on-going  and  proposed  research  at  the  Alberta  Research  Council,  the  University 
of  Alberta,  the  Network  of  Centers  of  Excellence,  and  forest  industries  throughout  the  Alberta  by 
evaluating  differences  among  dispersions  of  residual  materials  rather  thsm  evaluating  differences  among 
amounts  of  residual  materials.  All  harvest  patterns  evaluated  in  this  study  can  be  implemented  by  forest 
companies  that  leave  at  least  3%  of  the  canopy  trees  in  some  of  their  harvest  blocks.  The  project  started 
April  1997,  with  one  year  of  pre-harvest  biotic  sampling  and  will  continue  for  two  years  of  post-harvest 
sampling.  Although  the  primary  focus  of  the  project  is  on  the  biotic  community,  we  will  also  evaluate  the 
fall-down  rates  of  trees  and  snags  that  are  left  standing  in  the  experimental  areas. 

Within  both  study  areas  we  grouped  the  experimental  stands  into  four  stand  groups  to  increase  the 
statistical  power  of  the  tests.  Prior  to  harvest,  composition  and  density  of  tree  and  snag  species,  volume 
of  downed  woody  material,  bird  communities,  and  bat  abundance  varied  among  the  habitat  groups. 
However,  after  accounting  for  differences  among  habitat  groups,  forest  structure  and  biotic  communities 
did  not  V2iry  among  treatments.  Thus,  if  we  find  differences  after  harvest  they  can  be  attributed  directly 
to  the  experimental  treatments.  During  1998  and  1999  we  will  survey  the  forest  structure  and  biotic 
communities  in  these  experimental  stands  to  evaluate  the  relative  merits  of  leaving  big  versus  small 
patches  of  standing  trees  within  harvest  areas. 

Deliverables 

As  part  of  developing  and  conducting  this  research  we  have  worked  closely  with  government  regulatory 
agencies  and  the  forest  industry  to  ensure  that  the  results  will  be  applicable  to  these  organizations.  We 
will  continue  to  work  with  these  organizations  and  present  results  (to  groups  and  to  individuals)  as 
results  become  available.  Specific  deliverables  include: 

1)  Progress  reports  submitted  each  year. 

2)  Descriptions  of  how  the  spatial  pattern  of  standing  trees  and  snags  (i.e.,  variation  in  size  and  spacing 

of  tree  clumps)  affects  the  vegetation,  arthropod,  bird,  and  bat  communities  within  harvested  areas. 

3)  Descriptions  of  the  fall-down  rates  for  trees  and  snags  that  are  retained  in  harvested  areas.  Fall-down 

rates  will  be  evaluated  in  relation  to  the  size  and  placement  of  tree  clumps  in  harvested  areas. 

4)  Recommendations  for  clump  size  and  clump  spacing  that  are  most  suitable  for  meeting  biodiversity 

goals. 

5)  Presentations  of  research  findings  in  client  reports,  primary  journals,  and  at  conferences. 

6)  Participation  in  workshops,  planning  seminars,  and  committees  that  assess  and  regulate  alternative 

harvest  practices  in  Alberta. 
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CHAPTER  1.  BACKGROUND  AND  RATIONALE 


Background  and  Rationale 


Jim  Schieck 


Background 

As  a  consequence  of  having  a  thriving  forest 
industry  in  Alberta,  much  of  Alberta's  boreal 
forest  has  been  allocated  to  be  heirvested  either 
for  dimensional  lumber  or  for  pulp  and  paper. 
The  economic  benefits  of  this  intensive  forest 
harvesting  are  viewed  positively  by  Alberta's 
population.  The  ecological  ramifications  of 
extensive  forest  harvest,  however,  are  poorly 
understood  and  this  lack  of  knowledge  hampers 
management  decisions  and  land  use  planning. 
Consequently,  the  forest  industry  and 
government  organizations  in  Alberta  have  been 
conducting  research  to  develop  forest  practices 
that  are  ecologically  sound,  economically 
feasible,  and  sustainable. 

Short  rotation  practices  will  truncate  forest 
succession  over  much  of  the  landscape  and,  in 
some  areas,  may  result  in  old  forests  becoming 
relatively  scarce.  In  addition,  the  removal  of 
most  of  the  merchantable  trees  from  harvested 
areas  will  alter  natural  floral  and  faunal 
communities  in  those  areas  and  may  affect 
biotic  recovery.  Finally,  large  live  trees  and 
snags  provide  suitable  micro-habitats  for  many 
plants,  and  nesting  and  foraging  habitats  for 
many  animals,  but  with  conventional  harvest 
practices,  few  large  trees  and  snags  will  be 
present  in  harvested  landscapes. 

Management  strategies  are  being  developed  to 
produce  harvest  patterns  that  are  as  similar  as 
possible  to  natural  disturbances.  Since  large 
trees  and  snags  are  abundant  following  natural 
disturbances,  but  are  greatly  reduced  in 
harvested  areas,  forest  managers  have  been 
recommending  that  as  much  standing  residual 
material  as  is  socially  and  economically  feasible 
be  retained  within  harvested  areas.  Clumps  of 
standing  trees  and  snags  in  harvested  areas  may 
function  similar  to  skips  in  wildfire  areas,  and 
thus  act  as  refugia  for  some  species  while  the 
harvested  forest  re-grows.  If  residual  standing 
live  trees  continue  to  grow,  they  will  become 
large  live  trees  and  may  produce  some  large 
snags  in  mid-  and  late-rotation  forests.  These 


large  live  trees  and  snags  may  provide  some  of 
the  structures  and  micro-habitats  that  are 
normally  found  only  in  post-rotation  aged 
forests. 

Given  the  complexity  of  natural  ecosystems,  it 
is  not  surprising  that  the  ecological  and 
biodiversity  impacts  of  these  alternative 
harvest/silviculture  practices  are  poorly 
understood.  Only  by  experimenting  with  new 
techniques  in  the  field,  and  monitoring  the 
resulting  ecosystems  throughout  a  harvest 
rotation,  will  we  begin  to  understand  the 
potential  impacts  and  benefits  of  each 
technique.  To  this  end,  research  has  been 
conducted  by  scientists  within  the  Alberta 
Environmental  Protection,  the  University  of 
Alberta,  the  Network  of  Centers  of  Excellence, 
and  forest  industries  throughout  Alberta's  boreal 
forests.  Preliminary  results  from  past  and  on- 
going stand-level  research  indicates  that 
biodiversity  in  harvested  areas  is  positively 
related  to  the  amount  of  standing  trees  and 
snags  retained  in  those  areas.  However,  only 
anecdotal  information  is  available  to  evaluate 
whether  the  dispersion  of  those  standing 
residual  trees  also  affects  the  ecology  and 
biodiversity  in  harvested  areas.  The  present 
study  will  fill  that  gap  in  our  knowledge  by 
evaluating  whether  the  biota  present  in 
harvested  areas  with  small  clumps  of  trees 
differs  from  that  present  in  harvested  areas  with 
large  clumps  of  trees.  If  clump  size  has  a  direct 
impact  on  biota,  then  by  leaving  trees  and  snags 
in  the  appropriate  clump  sizes,  forest  managers 
may  obtain  benefits  to  biodiversity  at  little 
additional  cost.  The  proposed  experimental 
treatments  will  be  monitored  throughout  a 
harvest  rotation  (extended  monitoring  is  not  part 
of  the  present  proposal)  to  enable  researchers  to 
evaluate  changes  in  trees,  understory  vegetation 
and  wildlife  over  time.  This  controlled 
experimental  study  will  contribute  to  the 
adaptive  management  of  boreal  forests  and  will 
provide  information  that  forest  managers  can 
use  when  designing  harvest  plans. 
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Scientific  Rationale 

All  else  being  equal,  it  would  be  beneficial  to 
make  harvest  areas  as  similar  as  possible  to 
natural  disturbances  since  native  biota  evolved 
in  an  ecosystem  with  natural  disturbances 
(Urban  et  al.  1987).  The  dominant  natural 
disturbance  in  the  boreal  forests  of  North 
America  is  wildfires  (Eberhart  and  Woodard 
1987,  Hunter  1993),  and  cut-block  size  and 
shape  can  be  patterned  after  wildfires.  There  are 
limitations,  however,  to  matching  the  amount  of 
wood  remaining  in  post-harvest  areas  with  that 
remaining  post-wildfire  because  only  a  small 
percentage  of  the  wood  is  combusted  during  a 
wildfire  but  most  of  the  wood  is  removed  during 
harvest  (Lee  et  al.  1997).  Given  this 
discrepancy,  it  is  necessciry  to  manage  the 
relatively  few  trees  and  snags  that  are  retained 
in  harvest  areas  to  achieve  the  greatest  benefit 
for  biota.  Following  wildfires,  there  are  groups 
of  unbumed  trees  and  snags  (fire  skips) 
embedded  within  the  burned  areas  (Eberhart  and 
Woodard  1987).  These  fire  skips  could  be 
emulated  in  harvest  scenarios  by  leaving  clumps 
of  live  trees  and  snags  within  cut-blocks.  Both 
fire  skips  in  burned  areas,  and  residual  clumps 
of  trees  in  harvest  areas,  probably  retain  some  of 
the  structures  and  micro-habitats  that  were 
present  in  the  forest  prior  to  disturbance  and 
consequently,  may  act  as  refiigia  for  some 
native  species.  Small  clumps  may  be  sufficient 
for  a  few  species,  but  larger  clumps  probably  act 
as  refugia  for  more  species  (Saunders  et  al. 

1991)  .  For  example,  some  animals  require  large 
undisturbed  areas  for  protection  and  foraging 
(Pomeluzi  et  al.  1993),  and  some  plant  species 
require  large  undisturbed  areas  which  create 
suitable  micro-habitats  (Brothers  and  Spingam 

1992)  . 

The  refugia  potential  of  residual  clumps  of  trees 
is  probably  most  important  during  the  first  few 
decades  following  harvest  because  many  boreal 
forest  species  live  at  low  densities  in  young 
forests  (Schieck  et  al.  1995,  Kirk  et  al.  1996). 
Thus,  biota  that  are  present  in  the  residual 
clumps  may  disperse  from  these  into  harvest 
areas  once  the  new  cohort  of  trees  is  20  to  50 
years  old.  Distance  between  clumps  and  species 
dispersal  abilities  will  affect  the  rate  and 
dynamics  of  biotic  recolonization  in  harvested 


areas  (Saunders  et  al.  1991,  deMaynadier  and 
Hunter  1995).  The  degree  to  which  tree  clumps 
act  as  refugia  in  harvested  areas,  and  the 
appropriate  dispersion  of  those  clumps  to 
enhance  dispersal  of  biota,  have  not  been 
evaluated  in  boreal  forests. 

In  addition  to  the  refugia  potential  of  tree 
clumps,  there  may  be  management  value  to 
retaining  small  clumps  or  single  live  trees 
scattered  throughout  harvest  areas.  In  areas 
where  intensive  forest  harvest  will  decrease  the 
amount  of  old  forest  in  the  landscape,  it  may  be 
beneficial  to  produce  old-forest  structures  and 
micro-habitats  in  mid  serai  stages  of  forests 
(Spies  et  al.  1991).  Many  of  these  old-forest 
structures  and  micro-habitats  are  associated  with 
large  live  and  dead  trees  (Hansen  et  al.  1991, 
Schieck  et  al.  1995).  Consequently,  if  some 
merchantable  trees  are  retained  as  standing  live 
trees  at  harvest,  and  do  not  get  blown  down, 
they  will  become  large  by  mid-rotation  and  may 
accelerate  biotic  recovery  in  the  areas.  With  this 
management  strategy,  single  trees  or  many 
small  clumps  of  trees,  scattered  over  the  harvest 
area,  may  have  greater  impact  on  biota  than  if 
all  standing  trees  and  snags  had  been  retained  in 
few  large  clumps  because  old-forest  structures 
and  micro-habitats  would  then  be  dispersed 
throughout  the  harvest  area. 

The  most  appropriate  placement  of  standing 
trees  and  snags  differs  between  "managing  for 
refiigia"  and  "managing  for  the  production  of 
dispersed  large  trees".  Both  patterns  may  be 
useful,  but  only  by  experimentally  manipulating 
size  of  clumps  and  distance  between  clumps  and 
then  studying  the  biota  throughout  a  rotation 
will  we  understand  the  biological  and 
management  implications  of  each.  In  this  study, 
we  established  three  treatments  with  different 
clump  sizes  and  distances  between  clumps  in 
two  experimentally  harvested  areas  to  evaluate 
whether  biota  differ  among  the  treatments.  In 
the  present  study  we  will  evaluate  differences 
that  are  present  immediately  following  harvest. 
The  experimental  areas  will  be  maintained  for 
the  rest  of  the  harvest  rotation  however,  and 
proposals  to  evaluate  differences  among 
treatments  will  be  submitted  at  the  appropriate 
times  (i.e.,  once  the  new  cohorts  of  trees  are 
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well  established,  at  mid-rotation,  and  just  prior 
to  the  subsequent  harvest). 

Objectives 

In  the  present  study,  we  will  evaluate  how 
clump  size  and  distance  between  clumps  of 
residual  live  trees  and  snags  affects  biodiversity 
in  harvested  areas.  All  treatment  areas  will  be 
sampled  before  and  after  the  experimental 
harvest.  In  this  report  we  compared  pre-harvest 
biotic  communities  among  treatments  to  ensure 
that  treatments  were  similar  prior  to  harvest. 
Post-harvest  biotic  communities  will  be 
surveyed  during  1998  and  1999  and  compared 
among  treatments  to  evaluate  whether  the  biotic 
communities  present  within  the  first  two  years 
post-harvest  are  related  to  clump  size  and 
distance  between  clumps.  Our  results  will  be 
compared  to  that  found  in  benchmark 
communities  and  to  areas  where  other  harvest 
prescriptions  have  been  used. 

We  will  evaluate:  (1)  whether  the  bird 
community  that  has  specialized  foraging  and 
nesting  requirements  for  dead  trees  varies 
among  the  three  silviculture  treatments,  (2) 
whether  the  resident  bird  community  varies 
among  treatments,  (3)  whether  the  neotropical 
migrant  bird  community  varies  among 
treatments,  (4)  whether  the  bat  community 
varies  among  treatments,  (5)  whether  the 
arthropod  community  varies  among  treatments, 
(6)  whether  herb  and  shrub  communities  vary 
among  treatments,  and  (7)  whether  fall-down 
rates  of  standing  trees  and  snags  vary  among 
treatments. 
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CHAPTER  2.  STUDY  AREA  AND  METHODS 


Jim  Schieck 


The  study  is  being  conducted  in  two  forest  types 
in  the  boreal  forest  of  western  Alberta.  The 
study  area  southwest  of  Whitecourt,  AB  (54°  10' 
N,  116"  20'  W)  is  dominated  by  jack  pine  and 
black  spruce  that  is  approximately  1 10  years  old 
and  has  between  2,000  and  2,200  canopy  trees 
per  ha  (810-890  canopy  trees  per  acre). 
Throughout  the  area  there  is  localized 
topographical  relief  resulting  in  pockets  (<5% 
of  the  total  area)  of  aspen  and  white  birch. 
Twelve  stands  that  range  in  size  from  27  to  46 
ha,  mean  34  ha  (67  to  1 14  acres,  mean  84  acres) 
were  established  in  the  study  area  (Fig.  2.1). 
These  stands  were  categorized  into  four  groups 
of  three,  based  primarily  on  the  type  of 
vegetation  that  surrounded  them  and  secondarily 
on  the  tree  species  composition  in  the  area  prior 
to  harvest.  Within  each  stand  group  one  stand 
was  randomly  assigned  to  each  of  the  three 
treatments  (Fig.  2.1). 

Stands  were  harvested  between  September  1997 
and  February  1998.  Three  percent  of  the  pre- 
harvest  volume  was  left  in  tree  clumps  and  as 
many  of  the  non-merchantable  trees  and  snags 
as  was  possible  were  retained.  Clumps  were  not 
placed  on  roads  or  on  cut-lines,  but  otherwise 
they  were  dispersed  in  a  regular  pattern.  To 
reduce  wind-throw,  clumps  were  circular  and  all 
trees  within  3  m  (10  ft.)  of  the  perimeter  of  each 
clump  were  topped  at  4-5  m  (13-16  ft.)  to  create 
stubs.  We  have  assumed  that  each  stub  contains 
one  third  of  the  volume  of  the  original  tree  and 
have  added  sufficient  areaYo  each  clump  to 
replace  the  volume  that  was  removed  from  the 
stubs.  The  number  of  clumps  in  a  stand  was 
determined  based  on  the  treatment  being  created 
and  the  size  of  the  stand.  Based  on  a  stand  area 
of  30  ha  (74  acres),  treatment  one  would  have 
two  large  clumps  with  a  diameter  of  79.6  m 
(87.0  yd.)  and  each  clump  would  contain  850- 
935  residual  trees  and  145-160  stubs,  treatment 
two  would  have  eight  medium  sized  clumps 
with  a  diameter  of  41.7  m  (45.6  yd.)  and  each 
clump  would  contain  2(X)-220  residual  trees  and 
75-80  stubs,  and  treatment  three  would  have  32 
small  clumps  with  a  diameter  of  22.7  m  (24.8 


yd.)  and  each  clump  would  contain  45-50 
residual  trees  and  35-40  stubs.  During  the 
summer  of  1998  the  harvested  areas  within  the 
stands  will  be  re-planted  with  2-year  old  pine 
seedlings.  During  September  1999,  the  stands 
will  be  sprayed  with  herbicide. 

The  study  area  southwest  of  Grand  Prairie,  AB 
(54"  50'  N,  1 19"  20'  W)  is  dominated  by  aspen 
trees  that  originated  after  a  fire  approximately 
1(X)  years  ago.  This  smdy  area  has  a  few  small 
depressions  creating  pockets  of  white  spruce, 
balsam  poplar,  white  birch,  willow  and  black 
spruce.  Within  this  study  area  there  are 
approximately  800-1,000  canopy  trees  per  ha 
(324-405  Ccmopy  trees  per  acre).  The 
experimental  design  is  similar  to  that  for  the 
Whitecourt  study  area.  Twelve  experimental 
stands  were  created  by  harvesting  blocks  in  a 
checker-board  pattern  such  that  approximately 
half  of  the  smdy  area  Wcis  harvested  and  half  of 
the  area  unharvested  (Fig.  2.2). 

Stands  range  in  size  from  35  to  66  ha,  mean  52 
ha  (86  to  163  acres,  mean  128  acres).  These 
stands  were  categorized  into  four  groups  of 
three,  based  primarily  on  the  type  of  vegetation 
that  surrounded  them  and  secondarily  on  the 
tree  species  composition  in  the  area  prior  to 
harvest.  Within  each  group,  one  stand  was 
randomly  assigned  to  each  of  the  three 
treatments  (Fig.  2.2).  Stands  were  harvested 
between  September  1997  and  February  1998. 
Three  percent  of  the  pre-harvest  volume  was  left 
in  tree  clumps  and  as  many  of  the  non- 
merchantable  trees  and  snags  as  was  possible 
were  retained  in  each  stand.  Clumps  were  not 
placed  on  roads  or  cut-lines,  but  otherwise  were 
dispersed  in  a  regular  pattern  throughout  the 
stand.  The  number  of  clumps  in  a  stand  was 
determined  based  on  the  treatment  being  created 
and  the  size  of  the  stand.  Based  on  a  treatment 
area  of  50  ha  (124  acres),  treatment  one  would 
have  two  large  clumps  with  a 
diameter  of  97.7  m  (107.0  yd.)  and  each  clump 
would  contain  600  -  750  residual  trees, 
treatment  two  would  have  ten  medium  sized 
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FIGURE  2.1.  Locations  of  sites  and  stands  in  the  Whitecourt  study  area. 
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clumps  with  a  diameter  of  43.7  m  (47.8  yd.)  and 
each  clump  would  contain  120  -  150  residual 
trees,  and  treatment  three  would  have  50  small 
clumps  with  a  diameter  of  19.5  m  (21.3  yd.)  and 
each  clump  would  contain  24  -  30  residual  trees. 
No  post-harvest  silviculture  is  planned  for  the 
stands. 

To  ensure  that  harvest  patterns  were  consistent 
among  treatments,  the  boundaries  of  the  clumps 
of  trees  and  snags  that  were  retained  at  harvest 
were  marked  with  colored  plastic  ribbon  prior  to 
harvest.  No  machine  activity  occurred  within 
the  residual  clumps  of  trees  during  harvest. 

General  Survey  Methods 

Pre-harvest  surveys  of  biota  were  conducted 
during  the  spring  and  summer  of  1997.  Ten  1-ha 
(2.47  acre)  sites  were  established  in  the  center 
of  each  stand  (Figs.  2.1  and  2.2)  and  all  surveys 
of  forest  structure  and  biota  were  conducted  in 
those  central  sites.  Trees,  snags,  downed  woody 
material,  and  understory  vegetation  were 
surveyed  in  a  series  of  nested  plots  and  transects 
(Fig.  2.3,  see  Chapters  3  -  5).  Breeding  birds 
were  surveyed  using  a  combination  of  point 
counts  and  strip  transects  (Chapter  7).  Bats  were 
surveyed  using  automatic  remote  bat  detectors 
(Chapter  8).  Post-harvest  surveys  of  biota  will 
be  conducted  during  the  spring  and  summers  of 
1998  and  1999  (Table  2.1,  see  Chapters  3  -  8). 

Analysis  and  Reporting 

All  field  work  was  coordinated  by  Troy 
Sorensen  a  full-time  researcher  at  Alberta 
Research  Council.  Analyses  and  reporting  was 
conducted  during  the  fall  and  winter  by 
researchers  at  Alberta  Research  Council, 
Vegreville,  AB. 
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FIGURE  2.3.  An  example  of  the  transect  and  nested-plot  layout  for  a  1-ha  (2.47  acre)  site. 


TABLE  2.1.  Schedule  by  which  forest  structure  and  biota  will  be  surveyed  during  the  project 


Element 

1997 

1998 

1999 

(Pre-Harvest) 

(Post-harvest) 

(Post-harvest) 

Trees  and  Snags 

extensive  survey 

extensive  survey 

cursory  survey 

Down  woody  Material 

extensive  survey 

cursory  survey 

none 

Understory  Vegetation 

cursory  survey 

none 

extensive  survey 

Arthropods 

none 

extensive  survey 

cursory  survey 

Birds 

extensive  survey 

cursory  survey 

extensive  survey 

Bats 

extensive  survey 

extensive  survey 

extensive  survey 
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CHAPTER  3.  TREES  PRE-HARVEST  CHARACTERISTICS  AND  COMPOSITION  OF 
TREE  AND  SNAG  COMMUNITIES  IN  GRANDE  PRAIRIE  AND  WHITECOURT 
STUDY  AREAS 


Susan  Crites 


Introduction 

The  natural  disturbance  model  of  ecosystem 
management  assumes  that  biota  are  adapted  to 
natural  disturbances  and  are  more  likely  to 
maintain  viable  populations  if  forest  harvesting 
is  as  similar  as  possible  to  natural  disturbances 
(Swanson  and  Franklin  1992).  In  Alberta's 
boreal  mixedwood  forests,  fire  is  the  dominant 
natural  disturbance  and  is  currently  being  used 
as  a  template  for  timber  harvesting  by  some 
companies  (Anonymous  1997). 

Past  research  has  demonstrated  that  species 
respond  to  forest  structure  (Stelfox  1995). 
Timber  harvesting,  which  removes  structure, 
will  undoubtedly  have  negative  consequences 
on  biota  that  were  present  in  the  stand  before  it 
was  harvested.  Structure  in  the  form  of  large 
trees  and  snags  is  often  in  abundance  after  fire, 
and  its  presence,  mainly  in  the  form  of  snags 
and  downed  woody  material,  is  often  seen  long 
after  the  disturbance  (Lee  et  al.  1997).  To  allow 
timber  harvesting  to  more  closely  approximate 
fire,  forest  managers  are  advocating  that 
residual  merchantable  material  be  left  behind  at 
harvest.  The  appropriate  amount  and 
distribution  of  residual  material  on  cutblocks, 
however,  is  much  debated.  Clumps  left  behind 
at  harvest  could  be  compared  to  fire  skips  that 
serve  as  refugia  and  dispersal  sites;  therefore, 
larger  clumps  may  function  as  better  refugia. 
Conversely,  residual  material  evenly  distributed 
throughout  a  cutblock  could  serve  to  put 
structure  into  mature  forests  that  is  usually 
absent  from  harvest-origin  forests  until  older 
ages,  and  that  serves  as  habitat,  cover,  and 
foraging  sites.  Of  the  two  patterns,  it  is  not 
known  which  will  be  more  beneficial  to  biota. 

Trees  species  respond  differently  to  wind 
depending  on  their  structural  integrity,  decay 
properties,  exposure  to  wind  during 
development,  and  sheltering  ability  of 


surrounding  vegetation.  Trees  sheltered  within 
an  intact  forest  tend  to  have  smaller  boles, 
branches  concentrated  in  the  canopy,  and  are  not 
as  strongly  rooted  in  the  soil  as  trees  that  have 
been  exposed  to  wind  during  development. 
Newly  exposed  residual  trees  left  post-harvest 
are  susceptible  to  the  turbulent  and  accelerated 
wind  in  cutover  areas  because  they  have  not 
developed  the  structural  integrity  to  withstand 
such  winds  (Kimmins  1996). 

Deciduous  species  such  as  aspen  and  balsam 
poplar  are  more  likely  to  respond  to  wind  by 
snapping  off  along  the  bole  than  to  be  root 
thrown.  This  is  due  to  their  concentration  of 
branches  in  the  canopy,  the  depth  and 
complexity  of  their  rooting  system,  and  their 
susceptibility  to  fungal  attack  that  weakens 
boles  (Peterson  and  Peterson  1992).  Coniferous 
species  such  as  black  and  white  spruce  respond 
to  wind  by  being  root  thrown.  This  is  because 
of  their  shallow  rooting  system  and  high  surface 
area  (consisting  of  needles  and  branches) 
exposed  to  wind  that  extends  the  entire  length  of 
the  bole. 

Site  conditions  also  influence  the  ability  of  trees 
to  remain  standing  pre-  and  post-harvest.  In 
soils  saturated  with  water,  root  systems  are 
poorly  anchored  and  trees  have  the  tendency  to 
fall  over.  In  addition,  frost  heaving  can  uplift 
trees  from  the  soil.  Season  of  harvest  will 
influence  root  system  stability  and/or  tree  health 
in  areas  directly  affected  by  equipment.  Timber 
harvesting  when  the  ground  is  not  frozen  can 
lead  to  soil  compaction,  and,  if  compaction  is 
very  high,  it  can  kill  root  systems,  leading  to 
tree  death. 

In  this  report  I  evaluate  the  pre-harvest  density 
and  characteristics  (i.e.,  species,  heights,  and 
diameters)  of  trees  and  snags  among  stand 
groups  and  treatments  within  each  study  area 
(see  Chapter  2  for  experimental  design).  Post- 
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harvest,  I  will  determine  whether  tree  death  and 
tree  and  snag  fall  down  rates  vary  among 
treatments. 

Methods 

For  a  description  of  study  areas  and  methods, 
see  Chapter  2.  Trees  and  snags  were  Scmipled  in 
every  second  site  in  each  stand.  Each  one  ha 
sampling  site  was  divided  into  four  quadrants, 
and  alOxlOm(llxll  yd.)  tree  plot  was 
placed  in  two  randomly  selected  quadrants. 
Tree  plots  were  located  close  to  the  centre  of 
each  quadrant  (Fig.  2.3).  Within  each  plot, 
density,  diameter  at  breast  height  [1.4  m  (4.6 
ft.)],  height  class  [1  =  1.4-5  m  (4.6-16.4  ft.),  2  = 
5.1-10  m  (16.5-32.8  ft.),  3  =  10.1-20  m  (32.9- 
65.6  ft.),  4  =  >20  m  (>65.6  ft.)],  top  condition 
(broken  canopy,  broken  bole,  or  intact)  and 
species  were  recorded  for  trees  and  snags  that 
were  >1.4  m  (>4.6  ft.)  tall.  For  snags,  decay 
stage  was  also  recorded  (six  stages  from  Lee  et 
al.  1997;  stage  one  being  a  recently  killed  tree, 
stage  six  being  a  broken  off,  well-decayed  bole). 
Diameter  distributions  of  trees  and  snags  were 
plotted  and  evaluated  to  separate  stems  of  the 
canopy  cohort  from  stems  of  the  subcanopy  or 
understory.  Only  stems  of  the  canopy  cohort 
were  used  in  analyses.  For  snags,  "canopy" 
refers  to  trees  greater  than  a  particular  diameter 
(see  Results),  irrespective  of  whether  the  stem  is 
actually  part  of  the  canopy  (i.e.,  broken  off). 

General  Linear  Models  (GLM,  S  AS  Institute, 
Inc.  1989)  were  used  to  evaluate  differences  in 
abundances  of  canopy  trees  and  snags  among 
stand  groups  and  treatments.  Analyses  were 
completed  on  log  transformed  data.  Post  hoc 
tests  were  completed  using  Tukey's  Honestly 
Significant  Difference  test  (Zar  1984;  SAS 
Institute,  Inc.  1989).  All  results  were 
considered  statistically  significant  if  the 
probability  of  them  occurring  by  chance  was 
less  than  0.05. 

Results 

Live  Trees 

A  total  of  2672  and  4802  trees  and  snags  were 
sampled  in  Grande  Prairie  and  Whitecourt, 
respectively.  In  Grande  Prairie  there  were  four 
tree  species  found  during  the  surveys  with  aspen 
being  the  most  common  followed  by  balsam 


poplar,  paper  birch,  and  white  spruce  (Table 
3.1).  Seven  species  were  found  during  the 
surveys  in  Whitecourt,  with  jack  pine  and  black 
spruce  being  the  most  common  followed  by 
balsam  fir,  paper  birch,  aspen,  white  spruce,  and 
balsam  poplar  (Table  3.1). 

Stem  Diameter  and  Height 
Not  surprisingly,  stems  increased  in  diameter  as 
they  increased  in  height.  For  both  study  areas, 
the  majority  of  stems  were  short  [1.4-5  m  (4.6- 
16.4  ft.)]  and  in  diameter  class  0-5cm  (0-2  in.) 
(Figures  3.1  and  3.2).  In  Grande  Prairie,  the 
second  most  abundant  diameter  class  was  15.1- 
20  cm  (5.9-7.9  in.);  therefore,  canopy  trees  in 
Grande  Prairie  likely  consisted  of  trees  >15  cm 
(>5.9  in.)  dbh  (Figure  3.1).  The  majority  of 
trees  with  diameters  >15  cm  (>5.9  in.)  were  >20 
m  (>65.6  ft.)  tall  (Figure  3.1).  By  examining  the 
relationship  between  diameter  and  height  at 
Whitecourt,  Ccmopy  trees  likely  consisted  of 
trees  >10  cm  (>3.9  in.)  dbh  (Fig.  3.2).  Unlike 
trees  in  Grande  Prairie,  heights  of  canopy  trees 
were  between  10.1-20  m  (32.9-65.6  ft.),  with 
few  >20  m  (>65.6  ft.)  (Figure  3.2). 

Only  canopy  trees,  [i.e.,  stems  >10  cm  (>3.9  in.) 
and  >15  cm  (>5.9  in.)  dbh  for  Grande  Prairie 
and  Whitecourt,  respectively]  were  used  in 
subsequent  analyses.  Otherwise,  the  high 
frequencies  of  smaller  trees  in  both  study  areas 
would  dominate  density  calculations  and 
interpretations. 

Canopy  Stem  Density 

Average  number  of  canopy  stems  per  hectare 
was  746.8  (±75.0  SE)  [302.3  (±30.4  SE) 
stems/acre]  and  1441.8  (±123.5  SE)  [583.7 
(±50.0  SE)  stems/acre]  for  Grande  Prairie  and 
Whitecourt,  respectively.  In  Grande  Prairie, 
aspen  dominated  all  stand  groups  (Table  3.2). 
When  all  species  were  combined  in  Grande 
Prairie,  there  were  no  significant  differences  in 
stems  per  ha  among  stand  groups  or  treatments 
(Table  3.2).  For  individual  species  in  stand 
groups,  lower  densities  of  aspen  and  higher 
densities  of  balsam  poplar  were  found  in  stand 
group  C  (Table  3.2).  There  were  no  differences 
among  treatments  for  individual  species  (Table 
3.2). 
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TABLE  3.1.  Percent  species  composition  of  trees  and  snags  >1.4  m  (>4.6  ft.)  in  height  in 
Grande  Prairie  and  Whitecourt  study  areas.  Numbers  in  parentheses  indicate  the  actual  number 
sampled. 


Trees  Snags 


Species 

Grande  Prairie 

Whitecourt 

Grande  Prairie 

Whitecourt 

Populus  tremuloides 

81 (1937) 

2(79) 

78  (222) 

2(22) 

Populus  balsamifera 

8  (194) 

0(2) 

16(44) 

0 

Betula  papyrifera 

6(137) 

19  (726) 

4(12) 

9(87) 

Picea  glauca 

5(121) 

2(74) 

1  (2) 

1  (15) 

Pinus  banksiana 

0 

27(1033) 

0 

34  (335) 

Abies  balsamea 

0 

22  (828) 

0 

2(17) 

Picea  mariana 

0 

28(1070) 

0 

51 (507) 

unknown 

0 

0 

1(3) 

1(7) 

Total  (N) 

100  (2389) 

100(3812) 

100  (283) 

100  (990) 
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FIGURE  3.1.  Total  number  of  live  trees  categorized  by  diameter  and  height  classes  for  Grande 
Prairie. 
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FIGURE  3.2.  Total  number  of  live  trees  categorized  by  diameter  and  height  classes  for 
Whitecourt. 
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In  Whitecourt,  jack  pine  dominated  all  stand 
groups  and  treatments,  followed  closely  by 
black  spruce.  For  all  species  combined  in 
Whitecourt,  stand  group  A  had  the  highest  stem 
densities  and  stand  group  B  the  lowest  (Table 
3.3).  Jack  pine  and  black  spruce  were  present  in 
higher  densities  in  Stand  group  A,  and  in  lower 
densities  in  stand  group  B,  than  in  remaining 
stand  groups.  On  the  other  hand,  white  spruce 
and  balsam  fir  were  present  in  higher  densities 
in  stand  group  B  than  in  remaining  stand 
groups,  but  these  densities  were  relatively  low 
(Table  3.3).  When  all  stems  were  combined, 
there  were  no  significant  differences  among 
treatments  in  Whitecourt;  however,  aspen  was 
present  in  higher  densities  in  treatment  one  than 
in  remaining  treatments,  and  white  spruce  was 
absent  from  treatment  three  (Table  3.3). 

Snags 

Of  the  total  number  of  stems  sampled,  11% 
(283)  and  21%  (990)  were  snags  in  Grande 
Prairie  and  Whitecourt,  respectively.  Dominant 
snag  species  were  the  same  as  the  dominant  tree 
species  for  each  of  the  two  locations.  There 
were  no  paper  birch  snags  of  canopy  size  in 
Grande  Prairie,  and  no  balsam  poplar  snags  in 
Whitecourt. 

Snag  Diameter 

Diameter  distributions  for  snags  were  different 
than  for  live  trees  for  both  study  areas.  The 
snag  distribution  in  Grande  Prairie  peaked  at 
10.1-15  cm  (4-5.8  in.)  dbh  (Figure  3.3),  whereas 
in  Whitecourt,  the  distribution  was  high  from  0- 
10  cm  (0-3.9  in.)  dbh,  then  decreased  to  very 
few  trees  >30  cm  (>11.8  in.)  (Figure  3.4). 
Height  measurements  for  snags  are  not 
presented  because  many  snags  snapped  along 
the  bole.  The  same  diameter  cut  offs  as  used  for 
canopy  trees  was  used  for  the  snag  analyses 
(i.e.,  >15  cm  (>5.9  in.)  and  >10  cm  (>3.9  in.) 
dbh  for  Grande  Prairie  and  Whitecourt, 
respectively). 

Canopy  Snag  Density 

There  was  a  mean  of  77.5  (±19.2  SE)  [31.4 
(±7.8  SE)  per  acre]  and  289.3  (±44.8  SE)  [117 
(±18  SE)  per  acre]  large  snags  per  hectare  in 
Grande  Prairie  and  Whitecourt,  respectively.  In 
Grande  Prairie,  density  of  snags  did  not  v£iry 
significantly  among  groups  for  all  species 


combined  (Table  3.2);  however,  balsam  poplar 
and  birch  had  significantly  higher  densities  in 
stand  group  C  compared  to  remaining  groups 
(Table  3.2).  Density  of  large  snags  did  not  vary 
significantly  among  treatments  (Table  3.2). 

Similar  to  tree  density,  density  of  large  snags  in 
Whitecourt  varied  significantly  among  groups 
when  all  species  were  combined,  with  stand 
group  B  having  the  lowest  density  of  snags  and 
stand  group  A  the  highest  (Table  3.3),  The  low 
density  of  snags  in  stand  group  B  was  due  to  the 
low  densities  of  jack  pine  snags  (Table  3.3). 
Density  of  black  spruce  snags  was  significantly 
higher  in  stand  group  A  as  compared  to 
remaining  stand  groups  (Table  3.3).  Densities 
of  large  snags  did  not  vary  among  treatments  for 
all  species  combined,  but  for  individual  species, 
density  of  balsam  fir  snags  was  higher  in 
treatment  two,  than  in  remaining  treatments 
(Table  3.3). 

Decay 

In  all  dbh  classes,  snag  density  tended  to 
decrease  as  decay  state  increased  (Fig.  3.3  and 
3.4).  More  than  50%  of  the  snags  in  Grande 
Prairie  have  decay  stage  of  one,  with  82%  of  the 
snags  being  in  decay  stages  one,  two,  and  three 
(Figure  3.3).  In  Whitecourt  50%  of  the  snags 
were  within  decay  stages  one  or  two  and  95%  of 
snags  were  in  the  first  four  decay  stages  (Figure 
3.4). 

Top  Condition 

Of  all  snags  sampled,  66.2%  and  67.5%  had 
intact  boles  in  Grande  Prairie  and  Whitecourt, 
respectively.  In  Grande  Prairie,  26.6%  of  snags 
had  broken  boles  and  7.2%  had  broken 
canopies,  compared  to  16.6%  and  15.9%  in 
Whitecourt,  respectively. 

Discussion 

Composition  and  characteristics  of  the  tree  and 
snag  communities  were  different  between  the 
two  study  areas.  Grande  Prairie  was  an  aspen- 
dominated  community  with  canopy  trees 
commonly  reaching  heights  greater  than  20  m 
(65.6  ft.)  and  diameters  greater  than  15cm  (5.9 
in.)  dbh.  The  community  was  homogeneous 
with  only  four  species  and  having  aspen 
dominate  all  groups  and  treatments.  The 
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100 


Midpoint  of  Diameter  Class  (cm) 

FIGURE  3.3.  Toted  number  of  snags  categorized  by  diameter  class  and  decay  stage  for  Grande 
Prairie.  Decay  stage  one  refers  to  a  recently  dead  tree;  decay  stage  six  refers  to  a  well-decayed 
snag  that  is  soft  in  places  with  a  broken  bole. 


350 


ly/lidpoint  of  Diameter  Class  (cm) 

FIGURE  3.4.  Total  number  of  snags  categorized  by  diameter  class  and  decay  stage  for 
Whitecourt.  Decay  stage  one  refers  to  a  recently  dead  tree;  decay  stage  six  refers  to  a  well- 
decayed  snag  that  is  soft  in  places  with  a  broken  bole. 
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High  densities  of  trees  in  small  diameter  classes 
indicates  trees  are  continually  being  recruited 
into  the  understory  and  subcanopy.  If  gaps  form 
in  the  canopy,  these  understory  trees  may  grow 
to  become  part  of  the  canopy.  If  no  gaps  form 
before  these  understory  trees  are  out-competed, 
they  will  become  snags. 

As  with  trees,  aspen  also  dominated  the  snag 
community  in  Grande  Prairie  resulting  in  snags 
with  a  greater  frequency  of  broken  boles  than 
snags  in  Whitecourt.  The  self-thinning  nature 
of  aspen  stands  is  represented  in  the  median 
snag  diameter  of  10.1-15  cm  (4-5.9  in.)in 
Grande  Prairie,  which  is  one  diameter  cleiss 
lower  than  the  median  diameter  class  for  live 
canopy  trees  [15.1-18  cm  (6.0-7.1  in.)  dbh].  As 
aspen  stands  grow,  they  compete  for  resources 
and  continue  to  self-thin  until  they  reach  mature 
ages  and  canopy  equilibrium  is  obtained 
(Peterson  and  Peterson  1992).  As  a  result,  snags 
are  continuously  being  recruited  into  diameter 
classes  smaller  than  the  canopy  live  trees. 
Decay  stage  one  was  the  most  common  decay 
stage  for  all  snag  sizes,  with  the  exception  of 
snags  having  a  dbh  of  5.1-10  cm  (2.1-3.9  in.), 
and  snags  greater  than  30  cm  (1 1.8  in.)  dbh. 
Surprisingly,  many  snags  that  are  5.1-15  cm 
(2.1-5.9  in.)  dbh  have  been  present  long  enough 
to  have  decayed  to  stages  three,  four,  and  five. 

The  tree  community  in  Whitecourt  was  more 
heterogeneous  than  in  Grande  Prairie  with  two 
species,  jack  pine  and  black  spruce,  dominating 
and  eight  species  present.  Canopy  trees  in 
Whitecourt  were  shorter  and  had  smaller  dbh 
than  canopy  trees  in  Grande  Prairie.  As  in 
Grande  Prairie,  trees  in  Whitecourt  were 
abundant  in  small  diameter  classes,  indicating  a 
high  recruitment  of  trees  into  the  understory  and 
subcanopy.  The  fate  of  these  trees  will  depend 
on  canopy  gap  formation  and  competition 
dynamics  between  tree  species.  Tree  species 
richness  was  high  in  Whitecourt,  leading  to 
different  potential  successional  patterns  than 
found  in  Grande  Prairie.  Balsam  fir  and  paper 
birch  were  common  in  the  understory  of  many 
stands.  Balsam  fir  is  a  late  successional  shade- 
tolerant  species  that,  in  the  absence  of 
disturbance,  may  eventually  replace  the  current 
canopy  cohort  and  form  a  *climax'  community 
(Kimmins  1996). 


Based  on  vegetation  within  and  surrounding 
stands,  each  stand  was  assigned  to  one  of  four 
stand  groups.  Stand  groupings  accounted  for  a 
small  amount  of  variation  in  the  tree  and  snag 
communities  of  both  study  areas.  In  addition  to 
having  fewer  species,  the  tree  and  snag 
community  in  Grande  Prairie  exhibited  fewer 
differences  between  groups  than  in  Whitecourt, 
indicating  that  this  community  was  more 
homogeneous. 

There  were  no  differences  among  treatments  for 
tree  and  snag  communities  in  Grande  Prairie.  In 
Whitecourt,  there  were  differences  among 
treatments  for  aspen  and  white  spruce  live  trees, 
and  for  eispen  and  balsam  fir  snags.  When 
examining  the  Grande  Prairie  community  post- 
harvest,  any  differences  in  the  tree  and  snag 
communities  among  the  different  residual 
patches  will  be  attributed  to  the  harvesting 
treatment,  and  not  due  to  any  pre-harvest 
differences.  In  Whitecourt,  any  post-harvest 
differences  in  the  tree  and  snag  communities 
will  have  to  be  evaluated  after  taking  into 
account  the  differences  among  treatments  prior 
to  harvest. 

Future  Research 

Tree  and  snag  communities  will  be  surveyed 
and  monitored  in  1998  and  1999.  In  spring  of 
1998,  all  residual  trees  and  snags  will  be 
permanently  marked  within  a  subset  of 
randomly  selected  patches.  In  addition,  a  subset 
of  snags  and  trees  in  the  cut-over  area  will  be 
marked.  The  fate  of  this  residual  material  will 
be  regularly  determined  in  the  summer  of  1998 
and  1999.  In  1999, 1  will  determine  which 
configuration  of  patch  size  and  dispersion  is 
most  beneficial  to  the  long-term  persistence  of 
standing  trees  and  snags. 
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CHAPTER  4.  CHARACTERISTICS  OF  DOWNED  WOODY  MATERIAL  IN  j 
RELATION  TO  THE  DISPERSION  OF  RESIDUAL  TREES  FOLLOWING  HARVEST  | 

i 

Troy  C.  Sorensen 


Introduction 

Bryophytes,  lichens,  fungi,  invertebrates,  and 
herbaceous  plants  compose  the  majority  of 
biodiversity  in  an  ecosystem.  Many  of  these 
taxa  are  influenced  by,  or  dependent  on,  the 
abundance,  size,  spatial  distribution,  and 
decomposition  of  downed  woody  material 
(DWM;  McCullough  1948;  Thompson  1980; 
Harmon  et  al.  1986;  Soderstrom  1988a,b;  Crites 
and  Dale  1995,  in  press).  DWM  is  also  crucial 
as  nest  sites  (Thompson  1980),  food,  and 
protection  for  vertebrates  and  invertebrates 
(Harmon  et  al.  1986).  In  addition  to  its  direct 
influence  on  biodiversity,  DWM  functions  as  an 
important  reservoir  for  ecosystem  nutrients 
(Lambert  et  al.  1980,  MacMillan  1981,  Abbott 
and  Crossley  1982,  Harmon  et  al.  1986). 

Large  DWM  in  a  wide  rzuige  of  decay  stages  is 
often  stressed  as  important  in  maintaining  the 
species  composition  of  DWM-dependent  taxa 
(Soderstrom  1988a,b;  Lesica  et  al.  1991;  Crites 
and  Dale  1995,  in  press).  Clearcutting, 
however,  reduces  the  density,  volume,  size 
distributions  (Abbott  and  Crossley  1982;  Lesica 
et  al.  1991)  and  decay  class  distributions  of 
DWM  (Soderstrom  1988a,b;  Crites  and  Dale 
1995,  in  press;  Lee  et  al.  1995,  1997).  One 
strategy  to  minimize  these  effects  is  to  leave 
residual  live  trees,  snags,  and  DWM  in 
harvested  areas  (Lee  et  al.  1995,  1997)  in  an 
effort  to  promote  natural  size  and  decay  class 
distributions  of  DWM,  and  the  establishment  of 
DWM-dependent  taxa  (Soderstrom  1988b).  If 
forest  managers  are  willing  to  leave  this  residual 
material,  the  question  remains  of  how  to 
distribute  it  within  harvested  areas.  An  even 
spatial  distribution  of  DWM  will  facilitate  the 
dispersal  of  DWM-dependent  species,  but  will 
minimize  the  amount  of  interior-forest 
environment  which  promotes  natural 
decomposition  rates  of  DWM. 

Currently  this  study  has  two  objectives:  1)  to 
provide  pre-harvest  information  on  the  amount 


and  characteristics  DWM  in  the  two  study  areas 
and  2)  to  test  for  differences  in  DWM  among 
stand  groups  prior  to  harvest.  During  1998-99, 
post-harvest  density,  volume,  spatial 
distribution,  and  decomposition  of  DWM  will 
be  surveyed  and  compared  among  stand 
treatments.  Current  results  will  be  used  to 
account  for  pre-harvest  differences  among 
stands. 

Methods 

A  description  of  the  study  areas  and 
experimental  harvesting  is  presented  in  Study 
Area  and  Methods  (Chapter  2).  Ten  1-ha  (2.47 
acre)  sites  were  established  in  the  center  [ca, 
100  m  (109  yd.)  from  the  edge]  of  each  stand 
and  all  surveys  for  biota  were  conducted  within 
those  sites  (Fig.  2.1  and  2.2).  DWM  was 
surveyed  along  two  25  m  (27  yd.)  transects 
within  the  center  of  two  randomly  selected 
quadrants  of  each  site  (Fig.  2.3).  The  direction 
of  each  transect  was  random.  All  DWM 
intersecting  the  transect  with  diameters  >5  cm 
(>2  in.)  were  tallied;  in  addition,  the  diameter  at 
line  intersect  and  decay  class  (Table  4.1)  of 
large  pieces  [>10  cm  (>4  in.)]  were  measured. 
DWM  with  diameters  between  5-10  cm  (2-4  in.) 
and  >10  cm  (>4  in.)  will  be  called  small  and 
large  DWM,  respectively. 

The  12  stands  in  both  study  areas  were  assigned 
into  4  groups  to  group  stands  that  had  the  most 
similar  tree  composition.  The  three  stands  in 
each  group  were  then  randomly  assigned  a 
treatment  of  residual  tree  dispersion  (see 
Chapter  2).  Harvesting  occurred  during  the 
winter  of  1997-98. 

Linear  densities  of  DWM  were  calculated  as  the 
number  of  pieces  per  5  m  (16  ft.)  of  transect  and 
are  referred  to  in  the  rest  of  this  chapter  as 
densities.  The  volume  of  DWM  per  unit  area 
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TABLE  4.1.  Descriptions  of  decay  classes  (modified  from  Lee  et  al.  1997)  assigned  to  large 
[>10  cm  (>4  in.)  diameter]  DWM. 


Decay  Classes 

DWM  Characteristics 

1 

Log  whole  and  undecayed;  bark,  branches,  and  twigs  present  and 

llllCl^L,  iUg  CiCValCU  KjW  oUUL/Ull  UUllilo,  Yi\j\j\X  llOlVl. 

L 

Log  sound,  wood  hard;  twigs  mostly  lacking;  <50%  of  the  bark 

i 

Wood  soft  in  places;  some  branches  remaining;  >50%  bark  missing; 

wood  soft  in  places 

4 

Little  to  no  bark  remaining;  no  branches;  wood  fungi  present;  wood 

soft  with  small  crevices  and  small  pieces  lost. 

5 

Large  wood  fragments  lost;  outline  of  trunk  slightly  deformed; 

vascular  and  non- vascular  plants  beginning  to  colonize. 

6 

Wood  mostly  well-decayed;  log  colonized  by  various  vascular  and 

non- vascular  plants. 

7 

Humification  nearly  100%;  hard  to  define  as  log,  outline 

indeterminable;  no  evidence  of  hard  wood. 
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(v)  was  calculated  from  the  following  formula 
(Van  Wagner  1968): 

  1=1  


where  /  is  the  transect  length  [25  m  (27  yd.)],  d 
is  the  diameter  (cm)  of  the  large  pieces  [>10  cm 
(>4  in.)]  or  7.5  cm  (3  in.)  for  small  pieces  [5-10 
cm  (2-4  in.)],  and  /  is  the  number  of  pieces 
intersecting  the  transect.  The  resulting  units  are 
cubic  meters  per  hectare.  Densities  and 
volumes  were  calculated  for  both  small  and 
large  pieces  and  then  summed  for  a  total  value. 
Volume  was  also  calculated  for  each  decay  class 
(1-7)  of  the  large  pieces.  The  species 
composition  of  DW^  was  not  analyzed  because 
it  should  be  directly  related  to  the  species 
composition  of  trees  and  snags  (Chapter  3). 

General  Linear  Models  (GLM,  SPSS  1997) 
were  used  to  examine  differences  in  the 
densities  and  volume  of  DWM  between  stand 
groups  and  treatments.  Tukey's  Honestly 
Significant  Difference  tests  (a  =  0.05)  were 
used  for  all  post  hoc  tests  (Zar  1984;  SPSS 
1997).  Volumes  of  small  DWM  were  not 
analyzed  separately  because  they  were 
proportionally  related  to  small  DW^  densities. 
The  spatial  heterogeneity  of  DWM  distribution 
was  examined  by  using  a  nested  design  (sites 
within  stands). 

Results 

Grande  Prairie  Study  Area 

DWM  Densities 

There  was  0.99  ±  0.08SE  pieces  of  small  DWM 
and  0.94  ±  0.06  pieces  of  large  DWM  every  5m 
(16  ft.)  in  the  Grande  Prairie  study  area.  DW^ 
density  differed  among  stand  groups  for  small 
(F3  ,^  =  4.558,  p  =  0.005),  large  {F,  ,^  =  2.690,  p 
=  0.050),  and  total  DWM  (F,  ,^  =  6. 177,  p  = 
0.001 ;  Table  4.2).  Higher  densities  of  DWM 
were  found  in  stand  groups  near  bog  areas 
(Group  A  and  D,  Table  4.2,  Fig.  2.2). 

There  was  no  significant  difference  in  the 
densities  of  large  (F^^,^  =  0.027,  p  =  0.973)  or 


total  DWM  (F^,o8  =  1-930,  p  =  0.150)  between 
treatments;  however,  there  was  a  significant 
difference  between  treatments  in  the  density  of 
small  DWM  (F^,^  =  3.800,  p  =  0.025,  Table 
4.2). 

The  spatial  distribution  of  DWM  density  was 
heterogeneous  within  and  among  stands.  Small 
and  large  DW^  densities  were  spatially 
heterogeneic  both  within  (F,jjg  ,^0  ^  1.641,  p  < 
0.004)  and  among  stands  {F^^^  >  2.357,  p  < 
0.035).  On  the  other  hand,  total  DWM  density 
was  significantly  different  within  stands  (F,Qg 
=  1.927,  p  =  0.0002)  but  not  among  stands  (F^  ,^8 
=  1.757, /7  =  0.1 15). 

DWM  Volume 

There  was  77.4  m'/ha  (±  5.1SE)  [41.0  ydVacre 
(±  2.7SE)]  of  DWM  in  the  Grande  Prairie  study 
area.  Volume  differed  among  stand  groups 
(F3  ,^  =  2.952,  p  =  0.036,  Table  4.3)  with  group 
B  having  lower  total  DWM  volume  than  other 
groups  (Table  4.3);  potentially  due  to  their 
farther  distance  from  bogs  and  wet  areas  (Fig. 

2.2)  .  The  volume  of  large  DWM  did  not  differ 
among  stand  groups  (F^^^  =  2.487,  p  =  0.064, 
Table  4.3).  The  volume  of  large  and  total 
DWM  did  not  differ  among  stand  treatments 
prior  to  harvest  (F^,^  <  0.903,  p  >  0.408,  Table 

4.3)  . 

The  spatial  distribution  of  large  DWM  volume 
was  heterogeneous  both  within  (F^^^^  =1.531, 
p  =  0.012)  and  among  (F,  ,^  =  2.382,  p  =  0.034) 
stands;  however,  total  DWM  volume  was 
heterogeneous  within  stands  (F^^  ^^  =  1.584,/;  = 
0.007)  but  not  among  stands  F^  .^g'  =  1.981,  p  = 
0.075). 

DWM  Size 

The  average  diameter  of  large  DWM  in  the 
Grande  Prairie  study  area  was  15.5  cm  (± 
0.47SE)  [6.1  in.  (±0.19SE)].  However  the 
sizes  of  DWM  had  a  negative  exponential 
distribution  (Fig  4. 1). 

The  diameter  of  large  DW^  varied  among 
stand  groups  (F3  ,^  =  4.942,  p  =  0.003)  and 
treatments  (F^,„  =  3.415,  p  =  0.035;  Table  4.4) 
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TABLE  4.2.  Average  linear  density  of  DWM  (±  SE)  among  stand  groups  in  both  study  areas. 
Superscripts  denote  homogeneous  subsets  as  determined  by  Tukey's  Honestly  Significant 
Difference  (a  =  0.05)  for  significant  comparisons  among  groups/treatments. 


Linear  Density  (number  /  5m) 


Small  DWM        Large  DWM        Total  DWM 


Grande  Prairie  Area 

Group 

A 

1.24  ±0.18* 

1.09  ±0.10* 

2.33  ±  0.24  * 

B 

0  80  +  0  1 5  " 

0  78  4-  0  17  " 

1  58  +  0  10  " 

L.JO  ^  W.XW 

C 

0.90  ±0.06' 

0.89  ±0.11*' 

1.79  ±0.13'' 

D 

1.01  ±0.17*" 

0.99  ±0.11*' 

2.00  ±0.19*' 

Treatment 

1 

1.06  ±0.16* 

0.95  ±  0.02 

2.01  ±0.15 

2 

0.82  ±0.09' 

0.93  ±0.06 

1.75  ±0.15 

3 

1.09  ±0.14' 

0.93  ±0.18 

2.02  ±0.28 

Whitecourt  Area 

Group 

A 

1.16  ±0.09* 

1.16  ±0.06* 

2.32  ±0.11  * 

B 

0.85  ±0.19' 

2.02  ±  0.53  ' 

2.87  ±  0.63  ' 

C 

1.35  ±0.17* 

1.42  ±  0.20  *' 

2.77  ±  0.03  ' 

D 

1.08  ±0.09*' 

1.70  ±0.36*^ 

2.77  ±  0.37 ' 

Treatment 

1 

1.18±0.17 

L38±0.32* 

2.56  ±0.29 

2 

1.07  ±0.04 

1.74  ±0.43' 

2.82  ±0.45 

3 

1.08  ±0.20 

1.59  ±0.16*' 

2.67  ±0.11 
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TABLE  4.3.  Average  volume  per  unit  area  of  DWM  (±  SE)  among  stand  groups  in  both  study 
areas.  Superscripts  denote  homogeneous  subsets  as  determined  by  Tukey's  Honestly  Significant 
Difference  (a  =  0.05). 


Volume  (m  /  ha) 


Small  DWM 

Large  DWM 

Total  DWM 

Grande  Prairie  Area 

Group 

A 

17.3  ±2.5' 

62.7  ±  6.7 

80.0  ±4.7" 

B 

11.1  ±2.0' 

47.9  ±12.9 

59.0  ±10.8' 

C 

12.5  ±0.8' 

74.1  ±9.4 

86.7  ±9.2' 

D 

14.1  ±2.4^ 

69.8  ±11.3 

83.8  ±10.7' 

Treatment 

1 

14.7  ±2.2' 

63.5  ±  8.3 

78.2  ±7.5 

2 

11.3±1.3'' 

69.7  ±4.4 

81.0  ±5.3 

3 

15.2  ±2.0' 

57.7  ±  14.1 

72.9  ±13.8 

Whitecourt  Area 

Group 

A 

16.1  ±1.2' 

88.1  ±6.7' 

104.2  ±7.0' 

B 

11.8  ±2.6" 

239.6  ±  "8.4" 

251.4  ±58.9' 

C 

18.7  ±2.4' 

128.8-1-9.5' 

147.5  ±27.3' 

D 

14.9  ±1.2" 

164.3  ±54.0' 

179.2  ±  54.0' 

Treatment 

1 

16.3  ±  2.4 

136.8  ±43.1 

153.1  ±42.0 

2 

14.9  ±0.5 

170.8  ±  58.2 

185.7  ±58.6 

3 

14.9  ±  2.7 

158.0  ±35.7 

172.9  ±  33.2 
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with  larger  DWM  in  group  C  and  treatment  2 
(Table  4.4). 

DWM  Decay  Classes 

Very  little  of  the  large  DWM  was  undecayed 
(decay  class  1,  Fig.  4.2).  The  proportion  of 
DWM  lying  on  the  ground  (76.3%)  was  higher 
for  more  decomposed  DW^  (Fig.  4.2).  The 
decay  class  distribution  of  total  DW^  volume 
did  not  differ  among  stand  groups  (Fig.  4.3)  nor 
among  treatments  (Fig.  4.4). 

Whitecourt  Study  Area 

DWM  Densities 

There  was  1.11  ±  0.08SE  pieces  of  small  DWM 
and  1.57  ±  0.1 7SE  pieces  of  large  DWM  every 
5m  in  the  Whitecourt  study  area.  Density 
differed  among  stand  groups  (Table  4.2)  for 
small  (F3    =  4.233,  p  =  0.007),  large  (F,  = 
17.616,  p  =  2.22  X  10'),  and  total  DWM  (F,  ,^  = 
3.106,  p  =  0.030).  Low  densities  for  large  and 
total  DW^  occurred  in  stands  farthest  from 
riparian/creek  areas  (Group  A,  Table  4.2,  Fig. 
2.1). 

Varied  results  occurred  when  comparing  the 
densities  of  DW^  between  stand  treatments 
Table  4.2).  There  was  no  significant  difference 
in  the  density  of  small  {F^^^  =  0.466,  p  =  0.629) 
or  total  DWM  (F^^,^  =  1.138,  p  =  0.324)  among 
treatments;  however,  there  was  a  difference  in 
the  density  of  large  DW^  among  treatments 
(^1.08  =  5.646,  p  =  0.005,  Table  4.2).  These 
differences  are  probably  due  to  the  high  density 
of  large  DWM  in  B2  [2.97/5  m  (5.5  yd.) 
±0.1  BSE]  as  a  result  of  windfall  (pers.  obs.). 

Small,  large,  and  total  DWM  densities  were 
heterogenic  among  stands  (F^    >  2.445,  p  < 
0.030)  but  only  small  DW^  densities  were 
spatially  heterogeneic  within  stands  {F^^^^  = 
1.441,  p  =  0.026;  large  and  total  DWM,  p  > 
0.056). 

DWM  Volume 

There  was  170.6  m'/ha  (±  24.3SE)  [90.1 
ydVacre  (±  12.9SE)]  of  DWM  in  the  Whitecourt 
study  area.  DWM  volume  differed  among  stand 
groups  for  large  (F,  ,^  =  20.568,  p  =  1.29  x  10"'°) 
and  total  DWM  (F,;^  =  18.858,  p  =  6.58  x  10  '") 
with  group  A  having  the  lowest  volume  and 


group  B  having  the  highest  volume  (Table  4.3). 
Windfall  in  stand  B2  was  most  likely  the  cause 
for  the  extremely  high  volume  of  large  [330.7 
m'/ha  ±  20.7SE  (175.1  ydVacre  ±  1  LOSE)]  and 
total  DWM  [346.4  m'/ha  ±  20.8.SE  (183.4 
ydVacre  ±  1 1.0  SE)].  The  volume  of  large  (F^^j^g 
=  1.956,  p  =  0.146)  and  total  DWM  (F^,^  = 
1.759,  p  -  0.177)  did  not  differ  among 
treatments  (Table  4.3). 

The  spatial  distribution  of  large  and  total  DWM 
volume  was  heterogeneous  both  within  (F,og  > 
2.105,  p  <  4.01  X  10  ')  and  among  stands  (FJ,,^  > 
lL160,p<  1.16x10'). 

DWM  Size 

The  average  diameter  of  large  DW^  in  the 
Whitecourt  study  area  was  17.7  cm  (±  0.60SE) 
[7.0  in.  (±  0.24SE)];  however  the  sizes  of  DWM 
had  a  negative  exponential  distribution  (Fig. 
4.5).  Most  of  the  woody  material  (volume) 
constituted  pieces  with  diameters  of  10-30  cm 
(4-12  in.)  (Fig.  4.5). 

The  diameter  of  large  DWM  differed  among 
stand  groups  (F,  ,^  =  15.054,  p  =  1.418  x  IQ-*; 
Table  4.4)  with  group  A  having  the  smallest  and 
group  B  having  the  largest  DWM  (Table  4.4). 
Diameter  of  large  DWM  did  not  differ  among 
stand  treatments  (F^^^  =  0.170,  p=0.844). 

DWM  Decay  Classes 

Most  of  the  volume  of  DW^  was  in  advanced 
stages  of  decomposition  (decay  class  6  and  7, 
Fig.  4.6).  The  proportion  of  DWM  lying  on  the 
ground  (75.1%)  was  higher  for  more 
decomposed  DWM  (Fig.  4.6).  The  decay  class 
distribution  of  total  DW^  volume  did  not  differ 
among  stand  groups  (Fig.  4.7)  or  among 
treatments  (Fig.  4.8). 

Discussion 

Small  DWM  contributes  little  to  the  total  DWM 
volume  in  either  study  area;  however,  it  is 
preferred  over  larger  DWM  by  microarthropods 
(Abbott  and  Crossley  1982)  and  is  therefore  an 
important  factor  in  the  biodiversity  of  an 
ecosystem.  Furthermore,  small  DWM  is  an 
important  component  of  nutrient  cycling 
because  it  quickly  decomposes  (MacMillan 
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FIGURE  4.1.  Size  class  distribution  of  the  density  (histogram)  and  volume  (line)  of  DWM  in  the 
Grande  Prairie  study  area. 
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TAELE  4.4.  Average  diameter  of  large  DWM  (±  SE)  among  stand  groups  in  both  study  areas. 
Superscripts  denote  homogeneous  subsets  as  determined  by  Tukey's  Honestly  Significant 
Difference  (a  =  0.05). 


Diameter  (cm) 


Grande  Prairie  Area 

Liroup 

A 
/\ 

14.2  ±  0.5  / 

B 

15.011.39'"' 

C 

16.8  ±0.50' 

jj 

15.7  ±  0.26 

Treatment 

1 

15.510.80'"' 

2 

16.310.50* 

3 

14.611.03' 

Whitecourt  Area 

Group 

A 

15.510.2r 

B 

20.510.47" 

C 

17.210.63' 

D 

17.511.06' 

Treatment 

1 

17.91 1.09 

2 

17.41 1.04 

3 

17.711.31 
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FIGURE  4.2.  Average  linear  density  of  DWM  in  each  decay  class  in  the  Grande  Prairie  study 
area  distinguishing  between  DWM  raised  off  the  ground  and  DWM  lying  on  the  ground. 
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FIGURE  4.3.  Average  linear  density  (±  SE)  of  large  DWM  in  different  decay  stages  among  stand 
groups  in  the  Grande  Prairie  study  area. 
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FIGURE  4.4.  Average  linear  density  (±  SE)  of  large  DWM  in  different  decay  stages  among 
stand  treatments  in  the  Grande  Prairie  study  area. 


03/10/98 


Client  copy  -  not  for  citation  or  reproduction 


30 


Downed  Woody  Material 


FIGURE  4.5.  Size  class  distribution  of  the  density  (histogram)  and  volume  (line)  of  DWM  in  the 
Whitecourt  study  area. 
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FIGURE  4.6.  Decay  class  distribution  of  the  density  of  DWM  in  the  Whitecourt  study  area 
distinguishing  between  DWM  raised  off  the  ground  and  DWM  lying  on  the  ground. 
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FIGURE  4.7.  Average  linear  density  (±  SE)  of  large  DWM  in  different  decay  stages  among 
stand  groups  in  the  Whitecourt  study  area. 
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FIGURE  4.8.  Average  linear  density  (±  SE)  of  large  DWM  in  different  decay  classes  among 
stand  treatments  in  the  Whitecourt  study  area. 
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1981,  Abbott  and  Crossley  1982)  returning 
nutrients  to  the  soil. 

The  Grande  Prairie  study  area  had  higher 
density,  lower  volume,  and  smaller  sizes  of 
DWM  than  other  areas  of  >  120  years-old  aspen- 
dominated  forest  in  Alberta  (Lee  et  al.  1995, 
1997).  Densities,  volumes  and  sizes  of 
pine/spruce  DWM  could  not  be  found  in  the 
literature. 

The  clumped  spatial  distribution  of  DWM. 
within  stands  in  both  study  areas  suggests  that 
leaving  residual  material  in  patches  may  be 
closer  to  natural  spatial  distributions  (Lee  et  al. 
1995,  1997)  than  if  DWM  was  retained  in  a 
regular  pattern.  The  clumped  distribution  of 
DWTM  in  the  Whitecourt  study  area  is  probably 
due  to  the  nature  of  tree  mortality;  windfall 
usually  causes  clumped  distributions  of  DWM 
(Harmon  et  al.  1986).  In  both  study  areas, 
DW^  was  more  abundant  near  wet  areas 
possibly  because  the  roots  of  trees  in  wet  areas 
facilitate  windfall,  especially  for  coniferous 
species  (Harmon  et  al.  1986). 

The  decay  class  distribution  of  DWM  in  the 
Grande  Prairie  study  area  was  very  similar  to 
that  found  in  >  120  year-old  aspen  dominated 
forests  in  other  areas  of  Alberta  (Crites  and  Dale 
1995,  in  press;  Lee  et  al.  1995,  1997).  The 
onset  of  bark  loss  in  year  1  for  aspen  logs 
(Miller  1983)  results  in  DWM  passing  through 
decay  class  1  quickly.  Decay  class  distributions 
of  DWM  from  pine/spruce  dominated  forest 
could  not  be  found  in  the  literature. 

The  difference  in  decay  class  distributions  of 
DWM  between  the  Grande  Ptairie  and 
Whitecourt  study  areas  suggests  that  aspen, 
pine,  and  spruce  DWM  may  have  different 
residence  times  within  each  decay  class.  The 
shady,  cool  understory  of  the  pine/spruce  forest 
promotes  moss  cover  (30.61%  ±  1.28SE, 
Chapter  5)  compared  to  the  sunny,  dry 
understory  of  aspen  forest  (3.27%  ±  0.3  ISE, 
Chapter  5;  Soderstrom  1988a).  In  the 
pine/spruce  forest,  moss  quickly  covered  intact 
and  sound  DWM  therefore  most  DWM  was 
classified  as  decay  class  6  or  7.  On  the  other 
hand,  plants  may  colonized  DW^  at  similar 
times  in  both  smdy  areas  but  the  decomposition 


rate  of  DW^  following  plant  colonization 
(stages  6  -7)  may  be  slower  for  pine/spruce 
DWM  than  aspen  DWM.  The  time  required  for 
logs  to  lose  half  their  wood  density  varies 
greatly  between  soft-  and  hardwoods:  57  years 
for  pine  {Pinus  contorta,  Fahey  1983),  63  years 
for  balsam  fir  (Lambert  et  al.  1980),  but  only 
10-14  years  for  aspen  (Gosz  1080,  Miller  1983, 
reviewed  by  Harmon  et  al.  1986).  Hardwoods 
decompose  faster  than  softwoods  (MacMillan 
1981,  Harmon  1982)  due  to  different  wood 
structure  and  chemical  composition  (Wilcox 
1973). 

A  second  explanation  for  differences  in 
distributions  of  DWM  between  the  study  areas 
is  that  there  may  have  been  a  large, 
synchronized,  mortality  of  pine  and  spruce  trees 
in  the  WTiitecourt  study  area  within  the  time 
required  for  complete  decomposition  of  DW^. 
Balsam  fir  logs  can  retain  half  of  their  volume 
for  100  years  (Lambert  et  al.  1980).  This  cohort 
of  DWM  may  be  the  killed  trees  from  the  last 
fire  105  years  ago  (unpublished  data).  If  this  is 
true,  DWM  dynamics  not  reached  equilibrium 
before  the  area  was  harvested.  Within  a  few 
decades  the  DW^  in  decay  classes  6  and  7  may 
become  indistinguishable,  and  the  density  and 
size  distribution  of  DW^  in  the  Whitecourt 
study  area  may  be  simileu-  to  estimates  from  the 
Grande  Prairie  study  area  which  is  probably  at 
equilibrium. 

Natural  dynamics  in  old  forests  provide  constant 
inputs  of  DW^  resulting  in  a  relatively  even 
distribution  of  size  classes  and  decay  stages 
(Harmon  et  al.  1986;  Soderstrom  1988b;  Lee  et 
al.  1995,  1997;  Crites  and  Dale  1995,  in  press). 
On  the  other  hand,  managed  stands  have 
irregular  inputs  of  DW^  and  an  uneven 
distribution  of  DWM  decay  stages  (Soderstrom 
1988b).  Clearcutting  initially  results  in  lower 
density  and  volume  of  large  DWM  (Lesica  et  al. 
1991;  Lee  et  al.  1995,  1997)  and  a  pulse  of 
undecayed  small  branches  and  twigs  from  felled 
and  delimbed  trees  (Abbott  and  Crossley  1982; 
Crites  and  Dale  1995,  in  press).  At  mid-rotation 
the  small  residual  material  is  in  advanced  stages 
of  decomposition  with  low  input  of  trees  from 
self-thinning  and  death  of  residual  trees  (Lee  et 
al  1995,  1997).  In  old  age,  large  senescing  trees 
will  begin  to  consistently  input  DWM  with  even 
decay  class  and  size  distributions  (Lee  et  al. 
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1995,  1997).  Pulses  of  DWM  may  result  in 
temporal  heterogeneity  in  the  abundance  of 
species  dependent  on  specific  decay  stages  such 
as  epixylics  (Soderstrom  198^8a,b;  Lesica  et  al. 

1991)  and  of  species  dependent  on  specific  sizes 
of  DWM  such  as  microarthropods  (Abbott  and 
Crossley  1982,  Soderstrom  1988a). 

Many  species  of  bryophytes,  lichens,  and  some 
vascular  plants  are  dependent  on  specific  stages 
of  DW^  decomposition  (Muhle  and  LeBlanc 
1975;  Crites  and  Dale  1995,  in  press). 
Soderstrom  (1988b)  suggested  three 
requirements  to  maintain  pre-hcu^est 
compositions  of  species  dependent  on  DWM, 
such  as  lichens  and  bryophytes:  DWM  must  be 
present  in  all  stages  of  decomposition,  DWM 
must  have  an  even  spatial  distribution,  and  areas 
with  DWM  must  be  protected  from  drying  out. 
A  uniform  distribution  of  DWM  will  allow 
DW^-dependent  species  to  disperse  into 
harvested  areas  (Lesica  et  al.  1991).  Many  of 
the  rare  lichens  and  bryophytes  have  limited 
dispersal  capabilities  and  may  require  residual 
patches  of  DWM  to  act  as  dispersal  centers 
(Soderstrom  1988).  These  same  species  need 
shaded  areas  with  dependable  humidity  (Muhle 
and  LeBlanc  1975,  Larson  1984,  Soderstrom 
1988a,  Lesica  et  al.  1991,  Gustafsson  et  al. 

1992)  . 

The  smaller,  more  abundant  residual  patches 
created  in  treatment  1  of  this  study  will  result  in 
a  more  even  spatial  distribution  of  DW^  and 
may  facilitate  the  dispersal  of  DWM-dependent 
species  into  harvested  areas.  However,  the 
smaller  patches  may  have  more  variability  in 
evaporation  rates  which  suppresses  nonvascular 
plants  and  DWM  decomposition  (Gustafsson  et 
al.  1992,  Soderstrom  1988a).  Stands  in 
treatments  2  and  3  with  larger  and  wider  spaced 
residual  tree  patches  may  retain  more  interior- 
forest  habitat  with  a  more  natural  humidity, 
temperature,  shade,  and  natural  decomposition 
rates  of  DWM.  Natural  decomposition  rates 
may  maintain  temporally  homogenous 
compositions  of  DWM-dependent  species.  On 
the  other  hand,  the  dispersal  of  species  into  the 
harvested  areas  may  be  restricted  due  to  wider 
spaced  patches  of  DWM. 


Initially,  an  environment  that  allows  DWM- 
dependent  species  to  survive  will  be  more 
important  than  an  environment  that  facilitates 
dispersal.  The  ability  of  DWM-dependent 
species  to  disperse  will  become  important  once 
the  canopy  is  established  and  the  understory  of 
harvested  areas  are  buffered  from  sunlight  and 
desiccation.  Thus,  the  influence  of  DWM 
dispersion  on  biodiversity  may  shift  from  patch 
size  to  inter-patch  distance  as  the  harvested 
areas  regenerate. 

Pre-harvest  characteristics  of  DWM  provide  a 
benchmark  to  which  harvested  areas  can  be 
compared.  Stand  groups  were  successful  in 
accounting  for  pre-harvest  variance  of  DWM, 
whereas  there  were  few  differences  prior  to 
harvest  among  randomly  assigned  treatments. 
Thus,  any  post-harvest  differences  in  DWbA 
among  treatments  may  be  attributed  to  treatment 
effects. 

Future  Research 

Post-harvest  characteristics  of  DW^  will  be 
surveyed  in  1998  using  the  same  methods  as  in 
1997  with  additional  sampling  to  estimate 
DWM  characteristics  in  residual  tree  patches. 
Additional  transects  will  be  sampled  at  random 
directions  through  the  middle  of  residual  tree 
patches  for  a  total  of  80  m  (87.5  yd.)  per  stand. 

If  time  permits,  the  wood  density  of  large  DWM 
in  each  decay  class  will  be  sampled  so  that 
DWM  biomass  can  be  calculated  from  volume 
estimates.  Biomass  is  most  frequently  presented 
in  the  literature  (Gosz  1980,  Pastor  and 
Bockheim  1984,  Harmon  et  al.  1986).  For  each 
decay  class,  cross-sections  ('cookies')  of  ten  logs 
will  be  measured  for  volume  (diameter  and 
thickness  of  section),  weighed,  oven  dried,  and 
then  re-weighed.  This  will  also  provide 
estimates  of  moisture  content  and  biomass  loss 
which  in  turn  may  affect  bryophytes,  lichens 
and  other  plants  (Harmon  et  al.  1986). 

In  the  winter  of  1998-1999,  the  density,  volume, 
spatial  distribution,  and  decomposition  of  post- 
harvest  DWM  will  be  compared  among 
treatment. 
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CHAPTER  5.  A  SUMMARY  OF  UNDERSTORY  VEGETATION  IN  THE  BOREAL 
FOREST  PRIOR  TO  HARVEST. 


Christine  Gray 


introduction 

Due  to  funding  constraints,  only  preliminary 
surveys  were  conducted  for  the  understory 
vegetation  during  the  1997  field  season.  This 
report  is  a  summary  of  that  data.  More  detailed 
data  collection,  summary,  and  analyses  will  be 
conducted  during  1999/2000. 

l\/lethods 

In  both  study  areas,  vegetation  was  sampled  at  5 
out  of  the  10  possible  1  ha  (2.47  acre)  sites 
(every  second  site)  within  each  stand.  Each  1  ha 
site  was  divided  into  four  equal  quadrants  (Fig 
2.3).  Two  herb  plots  [0.5  x  0.5  m  (1.6  x  1.6  ft.)] 
were  systematically  placed  in  all  four  quadrants 
(8  herb  plots  per  site  and  40  herb  plots  per 
stand). 

For  sampling  purposes,  vegetation  was  divided 
into  two  height  strata.  The  herb  stratum 
included  all  herbs  regardless  of  height,  and  trees 
and  shrubs  that  were  <0.5  m  (<1.6  ft.)  in  height. 
The  shrub  stratum  included  all  shrubs  >0.5  m 
(>1.6  ft.)  in  height,  and  trees  that  were  between 
0.5  and  1.5  m  (1.6-4.9  ft.)  tall.  Shrubs  were 
divided  into  height  classes  of  0  =  0.5-1.39  m 
(1.6-4.6  ft.);  1  =  1.4-2.9  m  (4.7-9.5  ft.);  2  =  3.0- 
4.9  m  (9.6-16.1  ft.);  3  =  5.0-9.9  m  (16.2-32.5 
ft.);  4  =  10.0-19.9  m  (32.6-65.6  ft.);  and  5  = 
>20.0  m  (>  65.6  ft.).  In  the  Grande  Prairie  study 
area  one  shrub  plot  [2.0  x*2.0  m  (6.6  x  6.6  ft.)] 
was  systematically  placed  in  two  of  the  four 
randomly  selected  quadrants  (2  shrub  plots  per 
site  and  10  shrub  plots  per  stand).  In  the 
Whitecourt  study  area  a  shrub  plot  was  placed  in 
all  four  quadrants  (4  shrub  plots  per  site  and  20 
shrub  plots  per  stand).  Herb  plots  were  nested 
within  shrub  plots.  Sampling  was  conducted 
from  mid  June  through  mid  August  1997. 

Percent  of  ground  surface  cover  was  visually 
estimated  for  litter,  wood,  moss,  lichen,  and 
grass,  and  percent  of  ground  surface  covered  by 
a  vertical  projection  of  leaves  and  stems  was 


estimated  for  individual  species  of  herbs, 
shrubs,  and  trees.  In  the  Whitecourt  study  area 
percent  of  the  ground  surface  cover  by 
individual  moss  species  was  recorded  as  they 
were  the  dominant  ground  cover.  Grasses  and 
lichens  were  not  identified  to  species  at  either 
study  area.  Unknown  species  were  pressed  and 
later  identified  at  the  University  of  Alberta 
herbarium.  Nomenclature  follows  Moss  (1992) 
for  vascular  species.  For  nonvascular  species, 
nomenclature  follows  Vitt  et  al.  (1988). 

Results 

Grande  Prairie  Study  Area 

Eighty-three  species  were  found  in  the 
herbaceous  strata,  of  which  55  were  herbs,  25 
were  shrubs,  and  three  were  tree  species  (Table 
5.1).  Of  the  25  shrub  species  found  in  the  herb 
strata,  9  were  not  found  in  the  shrub  strata. 
Within  the  herb  strata,  much  of  the  ground  was 
covered  by  litter  (43%  cover)  although  herbs 
and  shrubs  were  also  abundant  (39%  and  18% 
cover  respectively).  Percent  cover  of  grasses, 
mosses,  and  lichen  were  much  lower  (Table 
5.1).  Bunchberry  was  the  most  abundant  herb 
(6%  cover),  with  Prickly  rose  and  Bracted 
honeysuckle  the  most  abundant  shrubs  (5% 
cover  for  each)  in  the  herb  strata.  Trees  had  low 
percent  cover  in  both  the  herb  and  shrub  strata 
(Table  5.1). 

The  shrub  strata  contained  21  species,  of  which 
18  were  shrubs  and  3  were  tree  species  (Table 
5.1).  Of  these.  River  alder  and  Common 
Labrador  tea  were  unique  to  the  shrub  stratum. 
The  most  abundant  shrubs  in  the  shrub  stratum 
were  Bracted  honeysuckle  (6%  cover), 
Canadian  buffaloberry  (5%  cover),  and  Prickly 
rose  (5%  cover).  Trembling  aspen  was  the  most 
abundant  tree  species  (1%  cover)  in  the  shrub 
stratum. 
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TABLE  5.1.  Mean  percent  cover  (±  SE)  of  understory  vegetation  for  the  Grande  Prairie 
study  area  during  1997.  Individual  species  are  grouped  by  strata  level  and  family.  


Species 

Mean 

Common  Name 

Species  Name 

Code 

Percent  Cover 

Litter  Total  Cover 

43.17  ±0.92 

Wood  Total  Cover 

5.31  ±0.45 

Lichen  Total  Cover 

0.01  ±  0.00 

Moss  Total  Cover 

3.27  ±  0.31 

Grass  Total  Cover 

4.81  ±0.26 

Herb  Total  Cover 

38.92  ±1.01 

Lycopodiaciaceae 

Stiff  club-moss 

Lycopodium  annotinum 

LYC  ANN 

0.06  ±  0.02 

Ground-cedar 

Lycopodium  complanatum  LYC  COM 

0.04  ±  0.04 

Equisetaceae 

Common  horsetail 

Equisetum  arvense 

EQU  ARV 

0.23  ±0.04 

Woodland  horsetail 

Equisetum  sylvaticum 

EQU  SYL 

0.04  ±0.01 

Polypodiaceae 

Lady  fern 

A  thyrium  filix-femina 

ATIi  HL 

0.01  ±0.01 

Oak  fern 

Gymnocarpium  dryopteris  GYM  DRY 

0.01  ±0.01 

unknown  fern 

UNK  FER 

0.08  ±0.08' 

Liliaceae 

Fairy  bells 

Disporum  trachycarpum 

DISTRA 

0.01  ±0.01 

Wild  lily-of-the-valley 

Maianthemum  canadense 

MAI  CAN 

0.67  ±0.08 

Star-flowered  false  solomon's-seal 

Smilacina  stellata 

SMI  STE 

0.04  ±  0.03 

Three-leaved  false  solomon's-seal 

Smilacina  trifolia 

SMITRI 

0.02  ±0.01 

unknown  lily 

Smilacina  spp. 

SMI  SPP 

0.02  ±  0.02 

Orchidaceae 

Rattlesnake-plantain 

Goodyera  oblongifolia 

GOO  OBL 

0.01  ±0.01 

Blunt-leaved  bog-orchid 

Habenaria  obtusata 

HAB  OBT 

0.01  ±0.01 

Ranunculaceae 

Red  and  white  baneberry 

Actaea  rubra 

ACT  RUB 

0.32  ±0.06 

Veiny  meadow  rue 

Thalictrum  venulosum 

THAVEN 

0.01  ±  0.01 

Saxifragaceae 

Bishop's  cap 

Mitella  nuda 

MITNUD 

0.78  ±0.08 

Rosaceae 

Wild  strawberry 

Fragaria  virginiana 

FRAVIR 

2.47  ±0.13 

Dewberry 

Rubus  pubescens 

RUB  PUB 

2.17  ±0.15 

Leguminosae 

Creamy  peavine 

Lathyrus  ochroleucus 

LATOCH 

3.01  ±0.18 

American  vetch 

Vicia  americana 

VICAME 

L20±0.11 

Violaceae 

Canada  violet 

Viola  canadensis 

VIOCAN 

0.26  ±0.07 

Kidney-leaved  violet 

Viola  renifolia 

VIOREN 

0.53  ±0.07 

Onagraceae 

Fireweed 

Epilobium  angustifolium 

EPIANG 

2.65  ±0.19 

'  Sum  of  two  unidentified  fern  species. 
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TABLE  5.1  cont. 


Species  Mean 

Common  Name  Species  Name  Code  Percent  Cover 

Araliaceae 


Wild  sarsaparilla 

Aralia  nudicaulis 

ARANUD 

1  83  ±0.25 

Umbelliferae 

Cow-parsnip 

Heracleum  lanatum 

HER  LAN 

0.23  ±0.09 

Comaceae 

Bunchberry 

Comus  canadensis 

COR  CAN 

6.11  ±0.32 

Pyrolaceae 

One-sided  wintergreen 

Orthilia  secunda 

ORT  SEC 

0.04  ±0.03 

Common  pink  wintergreen 

Pyrola  asarifolia 

PYR  ASA 

0.94  ±0.10 

Northern  starflower 

Trientalis  borealis 

TRIBOR 

0.02  ±0.01 

Gentianaceae 

Spurred  gentian 

Halenia  deflexa 

HAL  DEF 

0.06  ±  0.02 

Boraginaceae 

Tall  lungwort 

Mertensia  paniculata 

MER PAN 

1.61  ±0.17 

Labiatae 

Marsh  skullcap 

Scutellaria  galericulata 

SCU  GAL 

0.01  ±0.01 

Scrophulariaceae 

Common  red  paintbrush 

Castilleja  miniata 

CAS  MIN 

0.52  ±0.09 

Rubiaceae 

Northern  bedstraw 

Galium  boreale 

GAL  BOR 

1.36  ±0.09 

Sweet-scented  bedstraw 

Galium  triflorum 

GALTRI 

0.06  ±0.01 

Caprifoliaceae 

Twinflower 

Linnaea  borealis 

LIN  BOR 

2.89  ±0.17 

Compositae 

Common  yarrow 

Achillia  millefolium 

ACH  MIL 

0.26  ±  0.04 

Heart-leaved  arnica 

Amirn  rnrHif/'ilin 

ARNCOR 

0  36  +  0  07 

Lindley's  aster 

Aster  ciliolatus 

AST  CIL 

0.90  ±0.13 

Aster  conspicuus 

J.J  1  X  u.zy 

Purple-stemmed  aster 

Aster  puniceus 

AST  PUN 

0.03  ±0.02 

Narrow-leaved  hawkweed 

Hieracium  umbellatum 

HIE  UMB 

0.96  ±0.17 

Palmate-leaved  coltsfoot 

Petasites  palmatus 

PET  PAL 

2.38  ±0.16 

Arrow-leaved  coltsfoot 

Petasites  sagittatus 

PET  SAG 

0.14  ±0.08 

Canadian  butterweed 

Senecio  pauperculus 

SEN  PAU 

0.01  ±0.01 

Canada  goldenrod 

Solidago  canadensis 

SOL  CAN 

0.09  ±0.04 

Late  goldenrod  " 

Solidago  gigantea 

SOL  GIG 

0.03  ±  0.02 

Common  dandelion 

Taraxacum  officinale 

TAR  OFF 

0.01  ±0.01 

unknowns 

UNK 

0.09  ±0.07' 

Shrub  Total  Cover  (<  0.5  m) 

18.21  ±0.55 

Salicaceae 

willow 

Salix  spp. 

SAL  SPP 

0.10  ±0.04 

Betulaceae 

Green  alder 

Alnus  crispa 

ALNCRI 

0.03  ±0.02 

Beaked  hazelnut 

Corylus  comuta 

COR  COR 

0.07  ±0.05 

Grossulariaceae 

Skunk  currant 

Ribes  glandulosum 

RIBGLA 

0.01  ±0.01 

^  Sum  of  five  unidentified  herb  species. 
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TABLES.!  com. 


Species 

Mean 

Common  Name 

Species  Name 

Code 

Percent  Cover 

Northern  black  currant 

Ribes  hudsonicmum 

RIB  HUD 

0.06  ±  0.03 

Black  gooseberry 

Ribes  lacustre 

RIB  LAC 

0.04  ±  0.02 

Northern  gooseberry 

Ribes  oxyacanthoides 

RIBOXY 

0.1 6±  0.03 

Wild  red  currant 

Ribes  triste 

RIB  TRI 

0.36  ±  0.08 

unknown  currant 

Ribes  spp 

RIB  SPP 

0.10  ±0.04 

Rosaceae 

Saskatoon 

Amelanchier  alnifolia 

AME  ALN 

0.29  ±  0.08 

Pin  cherry 

Prunus  pensylvanicci 

PRU  PEN 

0.03  ±  0.02 

C^V\r\\c<^  f*li<»rrv 

PRTI VTR 

0  01  -J-  0  m 

Prir*Hv  roci* 
X  I  iwA^  y  iWow 

iXUdU  ULlLUlCiri^ 

RO<\  ACT 

J.J  J  -!-  U.^'T 

Wild  red  raspberry 

Rubus  idaeus 

u. /u  3:  u.  lu 

i^<llliJW~lCa.VCU  lllCO-UiJ WSWCCl 

opirucu  ciiua 

Ox  1  XJCd  1 

1  HQ  +  n  1 

Canada  hiiffalohpirv 

SHE  CAN 

0  62  +  0  09 

Comacpaf* 

Red-osier  dogwood 

{~' nmuv  vtnl/ini'Ffrn 

COR  STO 

0  12  +  0  05 

Ericaceae 

Common  bearberry 

Arrtn^tiinhvln^  iivfi-ur^i 

ARC  UVA 

0.01  ±0.01 

Common  labrador  tea 

Ledum  groenUindicum 

LED  GRO 

0.05  ±  0.03 

Velvet-leaved  blueberry 

Vaccinium  myrtilloides 

VAC  MYR 

0.78  ±0.14 

Bog  cranberry 

Vaccinium  vitis-idaea 

VAC  VIT 

0.04  ±0.02 

Canri  f ol  i  aceae 

TwinirnT  hr»np>\/«iipH<» 
X  vviiiiii^  iiuiiwy  dUwN.iw 

LAJfllLCrU  UlUlCU 

0  98  +  0  06 

iJlaL^lCU  llUIICjr  oUCKlC 

Lonicera  involucrata 

s  99  +  n  "^fi 
j,^^  ^  \j.j\j 

Common  snowberry 

Symphoricarpos  albus 

CVTVyf  AT  T5 

U.tKf  i  u.uz 

W  LIuoll-ClallUCliy 

Vibumunt  edule 

VTR  PDTT 

9  67  +  n  9Q 

Ay,jo  X  LmZfj 

willow 

Sn^lir  snr> 

SAL  SPP 

3.48  ±0.51 

Betulaceae 

vjiccii  aiuci 

Alnus  crispa 

ATM  PR  T 

1  M  -1-  0  69 

River  alder 

AlnuK  tpntiifnlin. 

ALN  TEN 

1.43  ±0.83 

Beaked  hazelnut 

/trvhiv  mmtitn 

COR  COR 

0.18  ±0.13 

Grossulariaceae 

Northern  gooseberry 

Ribes  oxyacanthoides 

RIB  OXY 

0.07  ±  0.04 

unknown  currant 

RIB  SPP 

0.02  ±  0.01 

Rosaceae 

Saskatoon 

Amelanchier  alnifolia 

AME  ALN 

0.94  ±  0.27 

Pin  cherry 

Prunus  pensylvanica 

PRU  PEN 

X  XX  \J  X  x««x  ^ 

0.41  ±0.16 

iink'nnwn  nhfrrx/ 

Prunus  spp. 

PRTT  ^PP 

Prunus  virginiana 

PRTT  VTR 
rlvU  V  lis. 

0  09  +  0  01 

Pnckly  rose 

Rosa  acicularis 

KUo  Ad 

^  no  -4-  n  ^ft 
j.yJZ  X  U.JO 

AVild  red  rasnhpirv 

RTfR  TDA 

0  14  +  0  07 

Vy«  1*T  .1.  \J»\J  i 

Elaegnaceae 

Canada  buffaloberry 

Shepherdia  canadensis 

SHE  CAN 

5.24  ±1.01 

Comaceae 

Red-osier  dogwood 

Comus  stolonifera 

COR  STO 

0.05  ±0.04 

Ericaceae 

Common  labrador  tea 

Ledum  groenlandicum 

LED  GRO 

0.08  ±0.08 
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TABLE  5.1  cont. 


Common  Name 


Species  Name 


Species 
Code 


Mean 

Percent  Cover 


Caprifoliaceae 

Twining  honeysuckle 
Bracted  honeysuckle 
Low  bush-cranberry 

Tree  Total  Cover  (<  0.5  m) 

Pinaceae 

White  spruce 
Salicaceae 

Trembling  aspen 
Betulaceae 

Paper  birch 

Tree  Total  Cover  (0.5  -  1.5  m) 

Pinaceae 

White  spruce 
Salicaceae 

Trembling  aspen 
Betulaceae 
  Paper  birch 


Lonicera  dioica 
Lonicera  involucrata 
Viburnum  edule 


Picea  glauca 
Populus  tremuloides 
Betula  papyrifera 

Picea  glauca 
Populus  tremuloides 
Betula  papyrifera 


LON  DIO 
LONINV 
VIB  EDU 


PIC  GLA 
POPTRE 
BET  PAP 

PIC  GLA 
POPTRE 
BET  PAP 


0.09  ±0.06 

6.27  ±1.14 

4.28  ±0.75 

0.07  ±  0.04 

0.01  ±0.01 
0.05  ±0.04 
0.01  ±0.01 
0.87  ±0.31 
0.04  ±0.04 
0.74  ±0.28 
0.08  ±0.08 
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Whitecourt  Study  Area 

Seventy-six  vascular  species  (including  53 
herbs,  19  shrubs,. and  4  tree  species)  and  22 
nonvascular  species  (19  mosses  and  3 
lichens)  were  found  in  the  herb  strata  (Table 
5.2).  Within  the  herb  strata  moss,  litter  and 
herbs  were  equally  abundant  (3 1  %,  30% 
and  28%  cover  respectively).  Knight's 
plume  (17%  cover)  and  Big  red  stem  (15% 
cover)  were  the  most  abundant  moss 
species.  Bunchberry  was  the  most  abundant 
herb  (4%  cover)  with  Common  Labrador  tea 
(3%  cover)  and  Velvet-leaved  blueberry 
(2%  cover)  the  most  abundant  shrubs  in  the 
herb  strata.  Balsam  fir  was  the  most 
abundant  tree  species  in  both  the  herb  and 
shrub  strata  (Table  5.2). 

The  shrub  strata  contained  a  total  of  20 
species,  of  which  14  were  shrubs  and  6  were 
trees  (Table  5.2).  Five  shrub  species  found 
in  the  herb  strata  were  not  found  in  the 
shrub  strata.  DeviFs  club  (2%  cover)  was 
the  most  abundcint  shrub  species  in  the 
shrub  strata.  Balsam  poplar  and  Trembling 
aspen  were  the  only  tree  species  found  in 
the  shrub  strata  and  not  the  herb  strata. 

Discussion 

The  Grande  Prairie  smdy  area  had  a  well 
developed  herb  and  shrub  strata;  shrubs 
were  tall  and  abundant  throughout  the  study 
area  and  the  herbs  covered  a  majority  of  the 
ground  space.  This  abundance  may  have 
been  due  to  the  deciduous  canopy  allowing 
light  to  penetrate  through  to  the  herb  and 
shrub  strata. 


The  Whitecourt  study  area  had  a  high 
density  of  coniferous  trees  (Chapter  4).  The 
herb  and  shrub  strata  were  less  developed  in 
this  area,  possibly  due  to  shading.  An 
abundant  moss  community  was  present, 
possibly  as  a  result  of  the  shaded  and  damp 
conditions. 

Due  to  funding  constraints  only  preliminary 
surveys  and  summary  were  conducted  for 
the  preharvest  understory  vegetation. 

Future  research 

Surveys  and  analysis  of  post  harvest 
understory  vegetation  will  take  place  during 
1999/2000.  Methods  will  be  similar  to 
those  presented  above  but  to  obtain 
adequate  data,  all  1-ha  sites  will  be 
surveyed  (ie.  twice  the  intensity  of  the 
surveys  during  1997).  Twenty  and  ten 
supplemental  herb  and  shrub  plots 
respectively  will  be  surveyed  within  each 
stand.  Half  of  the  supplemental  plots  will 
be  within  residual  tree  clumps  and  the  other 
half  >100  m  (>109  yd.)  from  residual  tree 
plots. 
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TABLE  5.2.  Mean  percent  cover  (±  SE)  of  understory  vegetation  for  the  Whitecourt  study  area 
during  1997.  Individual  species  are  grouped  by  strata  level  and  family.  


Species 

Mean 

Common  Name 

Species  Name 

Code 

Percent  Cover 

Litter  1  otal  Cover 

29.56  ±  1.08 

Wood  Total  Cover 

4.76  ±0.36 

Lichen  Total  Cover 

0.08  ±0.04 

club  lichen 

Cladonia  spp. 

CLA  SPP 

0.01  ±0.01 

Studded  leather  lichen 

Peltigera  aphthosa 

PEL  APH 

0.01  ±0.01 

Finger  felt  lichen 

Peligera  neopolydactyla 

PELNEO 

0.04  ±0.03 

Moss  Total  Cover 

30.61  ±  1.28 

sphagnum 

Sphagnum  spp. 

SPH  SPP 

3.74  ±0.68 

Common  hair-cap 

Polytrichum  commune 

POL  COM 

0.39  ±0.12 

Juniper  hair-cap 

Polytrichum  juniperinum 

T*/^T      TT  TXT 

POL  JUN 

0.30  ±  0.08 

Slender  hair-cap 

Polytrichum  strict um 

POL  STR 

0.49  ±0.1 2 

Common  northern  lantern  moss 

Cinclidium  stygium 

CIN  STY 

1.46  ±0.20 

Knight's  plume 

Ptilium  crista-castrensis 

PTICRI 

17.46  ±  1.03 

Big  red  stem 

Pleurozium  schreberi 

PLE  SCH 

14.75  ±0.91 

Stair-step  moss 

Hylocomium  splendens 

HYL  SPL 

0.80  ±0.14 

Broom  moss 

Dicranum  scoparium 

Die  SCO 

0.27  ±  0.06 

unknown  mosses 

MOSS 

0.56  ±  0.39' 

Grass  Total  Cover 

3.10  ±  0.36 

Herb  Total  Cover 

28.24  ±  0.95 

Lycopodiaciaceae 

Stiff  club-moss 

Lycopodium  annotinum 

LYC  ANN 

2.61  ±0.21 

Ground-cedar 

Lycopodium  complanatum 

LYC  COM 

0  03  +  0  02 

Equisetaceae 

Common  horsetail 

Equisetum  arvense 

EQU  ARV 

0.52  ±  0.08 

Woodland  horsetail 

Equisetum  sylvaticum 

EQU  SYL 

0.76  ±0.13 

Polypodiaceae 

Lady  fern 

Athyrium  filix-femina 

ATHFIL 

1.91  ±0.49 

Oak  fern 

Gymnocarpium  dryopteris 

GYM  DRY 

2.53  ±0.21 

Cyperaceae 

unknown  sedge 

Carex  spp. 

SED  SPP 

0.01  ±0.01 

Liliaceae 

Fairy  bells 

Disporum  trachycarpum 

DISTRA 

0.48  ±0.08 

Wild  lily-of-the- valley 

Maianthemum  canadense 

MAI  CAN 

0.22  ±0.04 

False  solomon's-seal 

Smilacina  racemosa 

SMI  RAC 

0.01  ±0.01 

Star-flowered  false  solomon's-seal 

Smilacina  stellata 

SMI  STE 

0.03  ±0.02 

Three-leaved  false  solomon's-seal 

Smilacina  trifolia 

SMITRI 

0.05  ±0.02 

unknown  lily 

Smilacina  spp. 

SMI  SPP 

0.02  ±0.01 

Indian  hellebore 

Veratrum  viride 

VER  VIR 

0.04  ±  0.03 

Orchidaceae 

Heart-leaved  twayblade 

Listera  cordata 

LIS  COR 

0.01  ±0.00 

'  Sum  of  ten  unidentified  moss  species. 
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TABLE  5.2  cont. 


Species 

Mean 

Common  Name 

Species  Name 

Code 

Percent  Cover 

Urticaceae 

Stinging  nettle 

Urtica  dioica 

URTDIO 

0.05  ±0.03 

Ranunculaceae 

Red  and  white  baneberry 

Actaea  rubra 

ACT  RUB 

0.05  ±0.03 

Goldthread 

Coptis  trifolia 

COP  TRI 

0.16  ±0.04 

Tall  larkspur 

Delphinium  glaucum 

DEL  GLA 

0.01  ±0.01 

Veiny  meadow  rue 

Thalictrum  venulosum 

THAVEN 

0.03  ±0.01 

Saxifragaceae 

Bishop's  cap 

Mitella  nuda 

MIT  NUD 

0.58  ±0.07 

Three-leaved  foamflower 

Tiarella  trifoliata 

TIATRI 

1.38  ±0.16 

Rosaceae 

Wild  strawberry 

Fragaria  virginiana 

FRA  VIR 

0.02  ±0.01 

Purple  avens 

Geum  rivale 

GEURIV 

0.02  ±0.01 

Cloudberry 

Rubus  chamaemorus 

RUB  CHA 

0.30  ±0.09 

Five-leaved  bramble 

Rubus  pedatus 

RUB  PED 

2.90  ±0.22 

Dewberry 

Rubus  pubescens 

RUB  PUB 

1.91  ±0.15 

Leguminosae 

Creamy  peavine 

Lathyrus  ochroleucus 

LAT  OCH 

0.01  ±0.01 

Balsaminaceae 

Spotted  touch-me-not 

Impatiens  capensis 

IMP  CAP 

0.01  ±0.01 

Violaceae 

Kidney-leaved  violet 

Viola  renifolia 

VIOREN 

0.19  ±0.04 

Onagraceae 

Small  enchanter's  nightshade 

Circaea  alpina 

CIR  ALP 

0.01  ±0.01 

Fireweed 

Epilobium  angustifolium 

EPI ANG 

0.93  ±0.11 

Araliaceae 

Wild  sarsaparilla 

Aralia  nudicaulis 

ARA  NUD 

2.36  ±0.25 

Umbelliferae 

Cow-parsnip 

Heracleum  lanatum 

HER  LAN 

0.06  ±0.06 

Comaceae 

Bunchberry 

Comus  canadensis 

COR  CAN 

4.34  ±0.20 

Pyrolaceae 

One-sided  wintergreen 

Orthilia  secunda 

ORT  SEC 

0.01  ±0.00 

Common  pink  wintergreen 

Pyrola  asarifolia 

PYRASA 

0.29  ±0.06 

Green  wintergreen 

Pyrola  virens 

PYR  VIR 

0.05  ±0.02 

Boraginaceae 

Tall  lungwort 

Mertensia  paniculata 

MERPAN 

0.29  ±0.07 

Scrophulariaceae 

Common  red  paintbrush 

Castilleja  miniata 

CAS  MIN 

0.01  ±0.01 

Rubiaceae 

Northern  bedstraw 

Galium  boreale 

GAL  BOR 

0.01  ±0.00 

Sweet-scented  bedstraw 

Galium  triflorum 

GAL  TRI 

0.06  ±0.02 

Caprifoliaceae 

Twinflower 

Linnaea  borealis 

LIN  BOR 

1.94  ±0.15 

Compositae 

Heart-leaved  arnica 

Arnica  cordifolia 

ARN  COR 

0.29  ±0.06 

Narrow-leaved  hawkweed 

Hieracium  umbellatum 

HIEUMB 

0.03  ±  0.03 

Palmate-leaved  coltsfoot 

Petasites  palmatus 

PET  PAL 

0.65  ±0.08 

Arrow-leaved  groundsel 

Senecio  triangularis 

SEN  TRI 

0.01  ±0.01 

Late  goldenrod 

Solidago  gigantea 

SOL  GIG 

0.03  ±0.03 

Common  dandelion 

Taraxacum  officinale 

TAR  OFF 

0.01  ±0.01 
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TABLE  5.2  cont. 


Species  Mean 

Common  Name  Species  Name  Code  Percent  Cover 


unknowns 

UNK 

0.06  ±  0.06' 

Shrub  Total  Cover  (<  0.5  m) 

10.63  ±  0.52 

Salicaceae 

willow 

Salix  spp. 

SAL  SPP 

0.01  ±0.01 

Betulaceae 

Green  alder 

Alnus  crispa 

ALN  CRI 

0.04  ±  0.02 

River  alder 

Alnus  tenuifolia 

ALNTEN 

0.01+0.01 

Grossulariaceae 

Northern  black  currant 

Ribes  hudsonianum 

RIB  HUD 

0.02  ±  0.02 

Northern  gooseberry 

Ribes  oxyacanthoides 

RIB  OXY 

0.47  ±  0.09 

Wild  red  currant 

Ribes  triste 

RIB  TRI 

0.11  ±0.04 

unknown  currant 

Ribes  spp. 

RIB  SPP 

0.06  ±  0.02 

Rosaceae 

Saskatoon 

Amelanchier  alnifolia 

AME  ALN 

0.09  ±  0.04 

Prickly  rose 

Rosa  asicularis 

ROS  ACI 

0.95  ±0.11 

Wild  red  raspberry 

Rubus  idaeus 

RUB  IDA 

0.33  ±  0.09 

Western  mountain  ash 

Sorbus  scopulina 

SOR  SCO 

0.0 1  ±0.01 

Araliaceae 

Devil's  club 

Oplopanax  horridum 

OPL  HOR 

0.91  ±0.26 

Ericeaceae 

Creeping  snowberry 

Gaultheria  hispidula 

GAU  HIS 

0.56  ±0.09 

Common  Labrador  tea 

Ledum  groenlandicum 

LED  GRO 

2.55  ±0.28 

Black  huckleberry 

Vaccinium  membranaceum 

VAC  MEM 

0.16  ±0.05 

Velvet-leaved  blueberry 

Vaccinium  myrtilloides 

VAC  MYR 

2.12  ±0.21 

Bog  cranberry 

Vaccinium  vitis-idaea 

VAC  vrr 

0.92  ±0.12 

Caprifoliaceae 

Bracted  honeysuckle 

Lonicera  involucrata 

LONINV 

0.60  ±0.1 

Low  bush-cranberry 

Viburnum  edule 

VIB  EDU 

0.72  ±0.11 

Shrub  Total  Cover  (>  0.5  m) 

7.00  ±0.72 

Salicaceae 

willow 

Salix  spp. 

SAL  SPP 

0.32  ±  0.09 

Betulaceae 

Green  alder 

Alnus  crispa 

ALN  CRI 

1.05  ±0.34 

River  alder 

Alnus  tenuifolia 

ALNTEN 

0.17±0.11 

Grossulariaceae 

Northern  black  currant 

Ribes  hudsonianum 

RIB  HUD 

0.04  ±  0.04 

Northern  gooseberry 

Ribes  oxyacanthoides 

RIB  OXY 

0.14  ±0.06 

Wild  red  currant 

Ribes  triste 

RIB  TRI 

0.01  ±0.01 

unknown  currant 

Ribes  spp. 

RIB  SPP 

0.03  ±0.02 

Rosaceae 

Saskatoon 

Amelanchier  alnifolia 

AME  ALN 

0.01  ±0.01 

Prickly  rose 

Rosa  acicularis 

ROS  ACI 

1.20  ±0.18 

Wild  red  raspberry 

Rubus  idaeus 

RUB  IDA 

0.16±0.13 

Araliaceae 

Devil's  club 

Oplopanax  horridum 

OPL  HOR 

1.57  ±0.44 

Ericaceae 

Common  Labrador  tea 

Ledum  groenlandicum 

LED  GRO 

1.06  ±0.24 

*  Sum  of  four  identified  herb  species. 
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TABLE  5.2  cont.  

Species  Mean 

 Common  Name  Species  Name  Code  Percent  Cover 


Caprifoliaceae 


Bracted  honeysuckle 

Lonicera  involucrata 

0  66  4-  0  1 Q 

Low  bush-cranberry 

Viburnum  edule 

VIB  EDU 

0.59  ±0.13 

1  ree  i  otai  cover  (,<  u.5  m ) 

A  zlO  4-  A  1  n 

Pinaceae 

ijaisani  iir 

Abies  bQisatneci 

ART  RAT 

\Vhite  spruce 

Piceo.  glaucQ 

n  01  +  O  01 

U.Ui  X  U.Ui 

Black  spruce 

Picea  mariana 

O  (Y)  -4-  n  0 1 
U.UZ  X  U.UI 

Betulaceae 

Pjir»pr  hirr*H 

Dciuni  fjujjyi  ijcf  u 

RFT  PAP 

0  07  +  0  0*^ 

iree  lotal  cover  (,U.5  —  1.5  m) 

1  Til  4-  A 
l./o  X  U.Z7 

i^indccde 

Balsam  fir 

Abies  balsamea 

ABI BAL 

1.03  ±0.23 

White  spruce 

Picea  glauca 

PIC  GLA 

0.01  ±0.01 

Black  spruce 

Picea  mariana 

PIC  MAR 

0.16  ±0.06 

Salicaceae 

Balsam  poplar 

Populus  balsamifera 

POP  BAL 

0.02  ±0.02 

Trembling  aspen 

Populus  tremuloides 

POPTRE 

0.03  ±0.02 

Betulaceae 

Paper  birch 

Betula  papyrifera 

BET  PAP 

0.53  ±0.13 
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CHAPTER  6.  ABUNDANCE  AND  SPECIES  COMPOSITION  OF  CARABID  BEETLES 
IN  RELATION  TO  THE  DISPERSION  OF  RESIDUAL  TREES  FOLLOWING 
HARVEST 


Troy  C.  Sorensen 


Introduction 

Although  arthropods  compose  a  large 
proportion  of  the  biodiversity  in  most 
ecosystems;  the  influence  of  forestry  practices 
on  these  taxa  is  often  overlooked.  Recent 
studies  have  suggested  that  ground  beetles 
(Coleoptera:  Carabidae)  may  be  sensitive  to 
forestry  practices  and  serve  as  ecological 
indicators  (Thiele  1977,  Butterfield  et  al.  1995). 
Arthropod  communities  are  influenced  by 
forestry  practices  because  they  are  sensitive  to 
temperature  (Ericson  1979),  moisture  (Niemela 
et  al.  1992,  Niemela  and  Spence  1994),  ground 
litter  (Greenslade  1964),  and  understory  and 
canopy  cover  (Niemela  and  Spence  1994). 
Furthermore,  the  dispersion  of  suitable  habitat 
patches  will  influence  community  composition 
because  dispersal  is  energetically  expensive 
(Atkins  1969)  and  species  differ  in  their  ability 
to  disperse  (Thiele  1977). 

Boreal  communities  of  carabid  beetles  have 
been  studied  in  relation  to  seasonal  changes 
(Niemela  et  al.  1992,  Digweed  et  al.  1995),  edge 
effects  (Halme  and  Niemela  1993,  Spence  et  al. 
1996),  fragmentation  (Halme  and  Niemela 
1993),  post-harvest  regeneration  (Spence  et  al. 
1996),  and  post-fire  regeneration  (Holliday 
1991,  1992;  Spence  et  al.  1-996).  The  main 
objective  of  this  project  is  to  determine  the 
influence  of  residual  tree  dispersion  on  the 
species  composition  and  abundance  of  carabid 
beetles  in  aspen-dominated  and  pine-spruce- 
dominated  boreal  forests  of  western  Alberta. 
Specifically,  I  will  examine  whether  the  carabid 
community  within  and  surrounding  residual 
patches  varies  among  treatments  (see  Chapter  2 
for  a  description  of  the  treatments),  and  (2) 
whether  differences  among  treatments  can  be 
explained  by  microhabitat  preferences  or 
dispersal  capabilities  of  specific  species. 


Future  Research 

Due  to  funding  constraints,  carabid  beetles  were 
not  surveyed  during  1997.  Surveys  from  four 
stands  of  unharvested  forest  within  each  study 
area  will  serve  as  benchmarks  for  post-harvest 
data. 

Carabid  beetles  will  be  surveyed  with  pitfall 
traps  because  the  traps  are  effective  (Obrtel 
1971,  Luff  1975,  Niemela  et  al.  1986,  Spence 
and  Niemela  1994,  Digweed  et  al.  1995),  and 
results  will  be  comparable  to  previous  studies 
(Spence  et  al.  1996).  Furthermore,  pitfall  traps 
sample  a  variety  of  other  taxa  including  ants  and 
spiders  that  could  be  identified  and  analyzed  at 
a  later  date  if  funding  becomes  available.  Pitfall 
traps  have  been  criticized  because  they  sample 
activity  rather  than  abundance  and  therefore 
catches  over  short  durations  are  dependent  on 
weather  and  the  period  of  life  cycle  (Ericson 
1979).  However,  if  traps  are  open  for  weeks  at 
a  time,  and  sampling  occurs  at  multiple  times 
during  the  snow-free  portion  of  a  year,  the  effect 
of  weather  will  be  negligible  and  catches  are 
correlated  with  density  (Baars  1979,  Ericson 
1979). 

Pitfall  traps  will  consist  of  a  removable  500ml 
(0.53  qt.)  plastic  container  set  inside  a  larger 
permsinent  container  [IL  (1.06  qt.)]  so  that 
repeated  sampling  will  not  disturb  the  soil  and 
litter  around  the  trap  (Spence  and  Nimela  1994, 
Digweed  et  al.  1995).  Each  trap  will  be 
partially  filled  2-3  cm  (0.8-1.2  in.)  with 
undiluted  ethylene  glycol  to  increase  trap 
effectiveness  and  preserve  the  specimens 
(Holopainen  1992).  A  15  x  15  cm  (5.9  x  5.9  in.) 
plywood  cover  will  be  suspended  with  nails 
over  each  trap  to  prevent  rain,  leaves,  and  litter 
from  falling  into  the  traps  (Spence  and  Niemela 
1994).  Trapping  will  occur  during  the  first  two 
weeks  of  each  month  from  early  May  to  the 
middle  of  September  for  a  total  of  5  trapping 
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sessions.  Traps  will  be  left  open  for  2  weeks  at 
a  time  with  a  2  week  recuperation  period 
between  trapping  sessions  (traps' closed)  to 
prevent  population  depletion. 

Traps  will  be  placed  along  radial  transects  from 
the  center  of  residual  patches  in  harvested 
stands,  and  along  similar  transects  in  the 
unharvested  stands  (Fig.  6.1).  Transects  will 
begin  at  the  center  of  the  patch  and  end  at  the 
half-way  point  between  patches  in  the  harvested 
areas.  The  number  of  traps  in  each  stand  will  be 
similar;  however  the  number  of  traps  used  to 
sample  each  patch,  and  the  number  of  patches 
sampled  within  each  stand  will  depend  on  the 
treatment  (Fig.  6.1).  Single  traps  will  be  placed 
no  less  than  20m  apart  to  insure  independent 
catches  (see  Digweed  et  al.  1995).  Previous 
studies  caught  fewer  (67%  less)  beetles  in  pine- 
dominated  than  aspen-dominated  forests 
(Spence  et  al.  1996);  to  compensate,  the  number 
of  traps  in  each  stand  will  be  higher  in  the 
Whitecourt  study  area  (24  traps)  than  the 
Grande  Prairie  area  (16-18  traps).  Carabid 
beetles  will  be  identified  to  species  according  to 
Lindroth  ( 1 96 1 - 1 969) .  Reference  collections 
are  available  at  the  Strickland  Museum, 
University  of  Alberta,  and  at  the  Northern 
Forestry  Centre,  Canadian  Forest  Service, 
Edmonton. 
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Figure  6.1.  Locations  of  pitfall  traps  for  the  three  patch  sizes  in  the  Whitecourt  Area.  Trapping 
design  is  similar  for  the  Grande  Prairie  study  area  with  the  following  distances  from  patch 
edges:  -50  m,  -20  m,  -10  m,  0  m,  20  m,  50  m,  120  m,  and  250  m  (-  55  yd.,  -  22  yd.,  -  1 1  yd.,  0 
yd.,  22  yd.,  55  yd.,  131  yd.,  and  273  yd.).  Traps  must  be  at  least  20  m  (22  yd.)  apart  to  insure 
statistical  independence.  Traps  will  be  placed  along  transects  in  unharvested  forest  with  the 
same  spacing  as  the  traps  in  large  patches. 
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CHAPTER  7.  VARIATION  IN  BIRD  COMMUNITIES  IN  RELATION  TO  THE 

DISPERSION  OF  STANDING  LIVE  TREES  AND  SNAGS  ' 

  \ 

-  '  -  ■  ,1 
Jim  Schieck 


Introduction 

As  a  result  of  government  regulations  (e.g. 
British  Columbia  Ministry  of  Forests  1995)  and 
exploratory  adaptive  management  (e.g.,  Alberta- 
Pacific  Forest  Industries  Inc.  1996)  many  forest 
companies  have  modified  harvesting  strategies 
by  leaving  standing  trees  and  snags  in  harvested 
areas.  By  providing  more  complex  habitat  than 
is  found  after  conventional  clear-cut  harvest, 
these  residual  trees  may  make  the  harvest  areas 
more  useable  by  native  birds  (Schieck  et  al. 
1995,  Norton  and  Hannon  1997).  In  addition,  by 
providing  inputs  of  large  snags  and  downed  logs 
throughout  succession,  residual  trees  may 
promote  recovery  of  biota  in  harvested  areas 
(Hobson  and  Schieck  submitted).  However, 
there  is  little  information  available  to  determine 
the  number  of  residual  trees  and  their  spacing 
that  will  maximize  recovery  of  native  birds.  As 
a  result,  forest  managers  are  subjectively 
deciding  those  patterns  (Alberta-Pacific  Forest 
Industries  Inc.  1996). 

One  possible  template  to  follow  when  leaving 
standing  trees  is  a  namral  dismrbance  model;  all 
else  being  equal,  making  harvest  areas  as 
similar  as  possible  to  natural  disturbances 
probably  will  be  beneficial  to  native  biota  since 
they  evolved  in  an  ecosystem  with  natural 
disturbances  (Urban  et  al.  1987).  Wildfire, 
disease,  and  wind-throw  events  all  create  natural 
disturbances  in  forests,  but  the  dominant  natural 
disturbance  in  the  boreal  forests  of  North 
America  is  wildfire  (Eberhart  and  Woodard 
1987,  Hunter  1993).  Thus,  wildfire  is  the  most 
realistic  template  on  which  to  pattern  harvest  in 
the  boreal  forest.  Following  wildfires,  there  are 
groups  of  unbumed  trees  (fire  skips)  embedded 
within  the  bumt-over  areas  (Eberhart  and 
Woodard  1987).  Fire  skips  may  retain  some  of 
the  structures  and  micro-habitats  that  were 
present,  and  may  act  as  refiigia  for  some  of  the 
birds  that  lived  in  the  forest  prior  to  the 
disturbance.  Fire  skips  could  be  emulated  in 
harvest  areas  by  leaving  clumps  of  live  trees 


within  cut-blocks.  Small  clumps  may  be 
sufficient  for  some  native  biota,  but  larger 
clumps  probably  act  as  refugia  for  more  species 
(Saunders  et  al.  1991,  MacKenzie  and  Steventon 
1996,  Seip  and  Parker  1997)  because  some 
animals  require  large  undisturbed  areas  for 
protection  and  foraging  (Pomeluzi  et  al.  1993, 
Schmiegelow  et  al.  1997). 

There  are  limitations,  however,  to  using  a 
wildfire  template  for  harvest  because  these 
disturbances  remove  different  amounts  of  wood 
from  the  system.  Only  a  small  percentage  of  the 
wood  is  combusted  during  wildfire  resulting  in 
many  snags,  and  when  these  snags  fall,  many 
pieces  of  downed  woody  material  are  present  in 
post-fire  areas  (Lee  et  al.  1997).  Following 
harvest,  however,  most  of  the  wood  is  removed 
from  the  area  and  the  trees  that  remain  are  live 
so  that  inputs  of  snags  and  down  woody 
material  is  lower  and  occurs  over  a  longer 
period  than  that  found  after  fire.  Given  this 
discrepancy,  managers  should  not  expect  to 
"mimic"  the  strucmre  retained  after  fires  and 
additional  management  strategies  may  be 
needed  (Hutto  1995).  It  may  be  useful,  on  at 
least  part  of  the  harvest  area,  to  retain  trees  in 
patterns  that  promote  specific  groups  of  native 
biota  that  would  otherwise  "fall  through  the 
cracks".  One  such  pattern  would  be  to  have 
single  live  trees,  or  small  clumps  of  live  trees, 
scattered  throughout  the  harvest  area.  If  these 
trees  do  not  get  blown  down,  they  will  become 
large  by  mid-rotation  with  the  result  that  many 
of  the  structures  and  micro-habitats  associated 
with  large  live  and  dead  trees  will  be  present  in 
mid-rotation  forests  rather  than  only  in  old 
forests  (Spies  et  al.  1991).  Large  live  and  dead 
trees  that  are  scattered,  although  uncommon 
following  wildfires,  may  accelerate  recovery  of 
native  biota  because  many  birds  use  such 
strucmres  (Hansen  et  al.  1991,  Schieck  et  al. 
1995).  With  this  management  strategy,  single 
trees  or  many  small  clumps  of  trees,  scattered 
throughout  the  harvest  area,  may  have  greater 
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positive  impact  on  native  biota  than  if  all 
standing  trees  and  snags  had  been  retained  as  a 
few  large  clumps. 

When  managing  for  refiigia,  standing  trees  are 
retained  in  a  different  pattern  than  they  are 
when  managing  for  the  production  of  dispersed 
large  trees.  Both  patterns  may  be  useful,  but 
only  by  experimentally  manipulating  size  of 
clumps  and  distance  between  clumps,  and 
studying  the  birds  living  in  those  areas  will  we 
understand  the  impacts  of  each.  In  this  study,  we 
establish  three  treatments  with  different  clump 
sizes  and  distances  between  clumps  in  two 
experimentally  harvested  areas  to  evaluate 
whether  bird  communities  differ  among  the 
treatments.  This  project  will  evaluate 
differences  that  are  present  immediately 
following  harvest.  The  experiment,  however, 
will  be  maintained  throughout  a  harvest  rotation 
and  the  value  of  each  the  three  treatments  will 
be  evaluated  at  various  times  in  the  future. 

Objectives 

I  compeired  pre-harvest  bird  communities 
among  stand  groups  and  treatments  to  evaluate 
whether  bird  communities  were  similar  in  the 
treatments  prior  to  harvest.  Specifically,  I 
evaluated:  (1)  whether  the  bird  community  that 
has  specialized  foraging  and  nesting 
requirements  for  large  trees  varied  among 
treatments,  and  (2)  whether  the  resident  and 
migratory  bird  communities  varied  among 
treatments.  This  is  the  first  year  of  a  3-year 
study.  Post-harvest  bird  communities  will  be 
compared  among  the  three  treatments  using  data 
that  will  be  collected  during  1998  and  1999.  In 
addition,  post-harvest  bird  communities  will  be 
compared  to  the  communities  found  in  other 
studies  within  naturally  disturbed  forests  and 
following  different  silviculture/harvest 
treatments. 

Methods 

A  description  of  the  study  areas  and 
experimental  harvesting  is  presented  in  the 
Study  Area  and  Methods  (Chapter  2). 
Descriptions  of  the  pre-harvest  trees,  snags, 
down  woody  materials,  and  understory 
vegetation  is  presented  in  Chapters  3-5. 


Pre-harvest  Bird  Surveys 
Depending  on  the  shape  of  the  stands,  either  two 
or  three  strip  transects  100  m  (109  yd.)  wide  and 
having  a  total  length  of  1  km  (0.62  mi.),  were 
established  in  each  stand  (Figs.  2.1  and  2.2). 
Strip  transects  were  divided  into  ten  1-ha  sites 
and  permanent  markers  were  located  at  the 
center  of  each  site.  To  avoid  censusing  birds 
that  lived  in  adjacent  habitats,  the  outer  edge  of 
the  strip  transects  were  a  minimum  of  70  m  (77 
yd.),  and  when  possible  >100m  (>109  yd.),  from 
adjacent  forest  types.  Birds  were  surveyed  in 
each  stand  during  four  12-day  periods  between 
9  May  and  26  June  1997.  During  each  survey, 
the  center  lines  of  the  transects  were  walked  at  a 
rate  of  100  m  per  10  min  between  first  light  and 
three  hours  past  sunrise  on  days  with  little  wind 
and  no  rain.  The  location  of  all  birds  heard  or 
seen  within  the  transects  were  recorded  on 
maps.  A  50  m  (54.7  yd.)  detection  distance  was 
used  because  this  represents  the  maximum 
distance  that  all  bird  vocalizations  could  be 
detected  reliably  (Emlen  and  DeJong  1981, 
Wolf  et  al.  1995,  Schieck  1997).  At  the  center  of 
each  1-ha  (2.47  acre)  site  a  5-minute  point  count 
was  conducted.  Movements  of  birds  were 
recorded  on  maps  and  summarized  so  that 
individual  birds  were  tallied  only  once.  When 
possible  unknown  birds  were  located  and 
identified  after  the  end  of  each  survey.  Bird 
vocalizations  that  could  not  be  identified  in  the 
field  were  recorded  using  cassette  tape  recorders 
and  parabolic  microphones  and  identified  later. 
The  order  that  stands  were  surveyed  during  each 
period  was  randomly  determined.  Observers 
were  trained,  and  methods  were  standardized 
prior  to  the  start  of  the  surveys. 

Analyses  of  Pre-harvest  Bird  Communities 

Bird  species  were  categorized  based  on  their 
migration,  nesting,  and  foraging  habits 
(Campbell  et  al.  1990,  Semenchuk  1992, 
Gauthier  and  Auby  1996).  Resident  and  migrant 
bird  species  may  be  affected  differently  by  the 
dispersion  of  residual  trees  because  residents 
live  in  boreal  forests  throughout  the  year.  In 
addition,  the  dispersion  of  standing  residual 
trees  may  have  a  greater  affect  on  birds  that  nest 
or  forage  in  large  live  or  dead  trees  than  those 
that  do  not.  Tree  Swallows  (see  Table  7.1  for 
scientific  names)  were  included  in  the  analyses 
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TABLE  7.1.  Abundance  (mean  ±  SE),  migration  habits,  and  nesting/foraging  habits  of  bird 
species  that  were  detected  in  the  pre-harvest  forest  in  the  Whitecourt  study  area  and  Grande 
Prairie  study  area  during'  1997. 


Common  Name  Genus  And  Species  Species    Migration       Use       Mean  Abundance 

Code       Habits      Big  Trees    (number/ 10  ha) 


Whitecourt  Study  Area 


Nothem  Goshawk 

Accipiter  gentilis 

NOGO 

Resident 

Yes 

0.17 

0.17 

Ruffed  Grouse 

Bonasa  umbellus 

RUGR 

Resident 

No 

0.08 

± 

0.08 

Common  Snipe 

Gallinago  gallinago 

COSN 

Migrant 

No 

0.17 

± 

0.11 

Yellow-bellied  Sapsucker 

Sphyrapicus  varius 

YBSA 

Migrant 

Yes 

0.67 

± 

0.31 

Downy  Woodpecker 

Picoides  pubescens 

DOWO 

Resident 

Yes 

0.08 

± 

0.08 

Hairy  Woodpecker 

Picoides  villosus 

HAWO 

Resident 

Yes 

0.33 

± 

0.14 

Three-toed  Woodpecker 

Picoides  tridactylus 

TTWO 

Resident 

Yes 

0.17 

± 

0.11 

Northern  Flicker 

Colaptes  auratus 

NOFL 

Migrant 

Yes 

0.50 

± 

0.26 

Alder  Flycatcher 

Empidonax  alnorum 

ALFL 

Migrant 

No 

0.08 

± 

0.08 

Least  Flycatcher 

Empidonax  minimus 

LEFL 

Migrant 

Yes 

0.08 

± 

0.08 

Tree  Swallow 

Tachycineta  bicolor 

TRSW 

Migrant 

Yes 

0.17 

± 

0.11 

Gray  Jay 

Perisoreus  canadensis 

GRJA 

Resident 

Yes 

1.58 

± 

0.5 

Black-capped  Chickadee 

Panis  atricapillus 

BCCH 

Resident 

No 

3.83 

± 

0.69 

Boreal  Chickadee 

Parus  hudsonicus 

BOCH 

Resident 

No 

2.92 

0.57 

Red-breasted  Nuthatch 

Sitta  canadensis 

RBNU 

Resident 

Yes 

1.17 

± 

0.71 

Brown  Creeper 

Certhia  americana 

BRCR 

Migrant 

Yes 

0.25 

± 

0.18 

Winter  Wren 

Troglodytes  troglodytes 

WIWR 

Migrant 

Yes 

0.08 

± 

0.08 

Golden-crowned  Kinglet 

Regulus  satrapa 

GCKI 

Migrant 

Yes 

3.50 

± 

0.62 

Ruby-crowned  Kinglet 

Regulus  calendula 

RCKI 

Migrant 

Yes 

7.25 

± 

0.58 

Swainson's  Thrush 

Catharus  ustulatus 

SWTH 

Migrant 

No 

2.75 

± 

0.41 

Hermit  Thrush 

Catharus  guttatus 

HETH 

Migrant 

No 

1.50 

± 

0.36 

American  Robin 

Tardus  migratoruis 

AMRO 

Migrant 

No 

0.17 

± 

0.11 

Cedar  Waxwing 

Bombycilla  cedrorum 

CEWA 

Migrant 

No 

0.17 

± 

0.11 

Solitary  Vireo 

Vireo  solitarius 

SOVI 

Migrant 

Yes 

1.58 

± 

0.43 

Warbling  Vireo 

Vireo  gilvus 

WAVI 

Migrant 

Yes 

0.58 

± 

0.19 

Red-eyed  Vireo 

Vireo  olivaceus 

REVI 

Migrant 

Yes 

0.33 

± 

0.22 

Orange-crowned  Warbler 

Vermivora  celata 

OCWA 

Migrant 

No 

0.25 

± 

0.13 

Yellow  Warbler 

Dendroica  petechia 

YEWA 

Migrant 

No 

1.08 

± 

0.34 

Yellow-nimped  Warbler 

Dendroica  coronata 

YRWA 

Migrant 

Yes 

10.58 

± 

0.54 

American  Redstart 

Setophaga  ruticilla 

AMRE 

Migrant 

No 

0.17 

± 

0.11 

Ovenbird 

Seiurus  aurocapillus 

OVEN 

Migrant 

No 

0.25 

± 

0.13 

Mourning  Warbler 

Oporomis  Philadelphia 

MOWA 

Migrant 

No 

0.33 

± 

0.14 

Common  Yellowthroat 

Geothlypis  trichas 

COYE 

Migrant 

No 

0.08 

± 

0.08 

Canada  Warbler 

Wilsonia  canadensis 

CAWA 

Migrant 

No 

0.25 

± 

0.25 

Rose-breasted  Grosbeak 

Pheucticus  ludovicianus 

RBGR 

Migrant 

No 

0.08 

± 

0.08 

Chipping  Sparrow 

Spizella  passerina 

CHSP 

Migrant 

Yes 

5.33 

± 

0.63 

Lincoln's  Sparrow 

Melospiza  lincolnii 

LISP 

Migrant 

No 

0.08 

± 

0.08 

White-throated  Sparrow 

Zonotrichia  albicollis 

WTSP 

Migrant 

No 

2.25 

± 

0.62 

Dark-eyed  Junco 

Junco  hyemalis 

DEJU 

Migrant 

No 

4.50 

± 

0.66 

Brown-headed  Cowbird 

Molothrus  ater 

BHCO 

Migrant 

No 

0.25 

± 

0.13 

White-winged  Crossbill 

Loxia  leucoptera 

WWCR 

Resident 

Yes 

0.08 

± 

0.08 

Pine  Siskin 

Carduelis  pinus 

PISI 

Migrant 

Yes 

1.08 

± 

0.36 
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TABLE  7.1.  cont. 


Common  Name  Genus  And  Species  Species    Migration       Use       Mean  Abundance 

Code       Habits      Big  Trees    (number/ 10  ha) 


Grande  Prairie  Study  Area 


Ruffed  Grouse 

BonosQ  utnbcllus 

KUvJK. 

tvesiueni 

INO 

O  14 
U.  I'f 

Common  Snipe 

Gallituigo  gcillincigo 

Migrant 

INO 

u.uo 

U.Uo 

Yellow-bellied  Sapsucker 

Sphyrapicus  varius 

VT5  C  A 

Migrant 

1  es 

J. ZD 

— 

n  /1 1 

U.'H 

Downy  ^Voodpecker 

Picoidcs  pubcsccns 

lJ\J  W  w 

i\.esiocni 

I  CS 

yj.L  / 

n  17 

U.  1  / 

Hairy  Woodpecker 

Picoides  villosus 

TT  A 

rLA.  W  VJ 

Resident 

I  es 

KJ.JJ 

± 

u.iy 

Northern  Flicker 

CoUiptcs  ourotus 

iNVJrL, 

Migrant 

I  es 

u.o  / 

U.jO 

r  iieaieci  w  ooopecKer 

uryocopus  piieutus 

PTwn 

ivcsiucni 

I  Cd 

0  OR 

it 

n  OR 

U.Uo 

Alder  Flycatcher 

EtnpidoTicix  olnorum 

AT  T7T 

Migrant 

Mo 

u.oj 

0  71 
U.Zl 

i-icobi  Fiycaicner 

Eitipidonox  tninimus 

jvngrani 

I  es 

R 

O.J  J 

1  04. 

Tree  Swallow 

Tachycineta  bicolor 

IKoW 

Migrant 

1  es 

U.l  / 

O  17 
U.l  / 

Gray  Jay 

Perisoreus  canadensis 

nx>  T  A 

Resident 

I  es 

— 

u.Zj 

Common  Raven 

Corvus  corax 

Resident 

I  es 

\J.ZD 

O  1  fi 
U.io 

Black-capped  Chickadee 

Panis  atricapillus 

Resident 

INO 

Z.Jo 

O 

U.JO 

Boreal  Chickadee 

Parus  hudsonicus 

Resident 

iNO 

U.Zj 

nil 

Red-breasted  Nuthatch 

Sitta  canadensis 

DT5MTT 
rvDlNU 

Resident 

1  es 

U.JO 

O  7Q 

u.zy 

Golden-crowned  Kinglet 

Regulus  satrapa 

Migrant 

Yes 

U.Uo 

U.Uo 

Ruby-crowned  ICinglet 

Regulus  calendula 

Migrant 

I  es 

n  AO 

± 

O  7A 

u.zo 

Veery 

Catharus  fuscescens 

VtiriK 

Migiant 

INO 

U.Uo 

— 

U.Uo 

Swainson's  Thrush 

Cathanis  ustulatus 

oW  In 

Migrant 

INO 

o.Uo 

— 

U.4j 

riermii  inrusn 

Catharus  guttatus 

rltlrl 

Migrant 

INO 

1  HQ 
l.Uo 

± 

n  7A 
U.ZO 

American  Robin 

Turdus  ntigratoruis 

AJvlKVJ 

Migrant 

INO 

l.O/ 

U.JJ 

Cedar  \Vaxwing 

Bombycilla  cedrorum 

Migrant 

INO 

U.ZD 

± 

n  1  fi 
U.lo 

Solitary  Vireo 

Vireo  solitarius 

oUVl 

Migrant 

V/»o 

Yes 

U.Zj 

— 

O  1  s 
U.lo 

Warbling  Vireo 

Vireo  gilvus 

WAVl 

Migrant 

Yes 

l.ZD 

± 

U.J  / 

Red-eyed  Vireo 

Vireo  olivaceus 

KbVl 

Migrant 

Yes 

J.yz 

— 

n  AO 
u.oy 

Tennessee  Warbler 

Vermivora  peregrina 

lew  A 

Migrant 

INO 

O  R7 

u.oz 

Orange-crowned  Warbler 

Vermivora  celata 

VJL,W  A 

Migrant 

INO 

n  S8 

U.JO 

U.ZJ 

Connecticut  Warbler 

Oporomis  agilis 

/~'r\\\T  A 
L.UWA 

Migrant 

INO 

l.Uo 

± 

O  70 

u.zy 

Yellow  Warbler 

Dendroica  petechia 

VT3W  A 
I  EWA. 

Migrant 

Mo 
INO 

R'^ 

0  74. 

U.  /** 

Magnolia  Warbler 

Dendroica  magnolia 

\yf  A  W  A 
JVIAW  A 

Migrant 

Mo 
INO 

n  ^R 

U.JO 

0  9"^ 

U.ZJ 

Yellow-rumped  Warbler 

Dendroica  coronata 

V1?\X/ A 
1  KW  A 

Migrant 

les 

Q  <in 
y.ju 

O  R7 
U.o  / 

Black-throated  Green  Warbler  Dendroica  virens 

r>  i  vjw 

Migrant 

I  es 

U.JJ 

0  99 
u.zz 

Blackpoll  Warbler 

Dendroica  striata 

DrWA 

Migrant 

INO 

n  OR 

U.Uo 

0  OR 

U.UO 

Black  and  white  Warbler 

Mniotilta  varia 

D  W  W  rV 

xviigrani 

Mo 

iNO 

0  OR 

0.08 

American  Redstart 

Setophaga  ruticilla 

/\iVllvC 

Migrant 

Mo 
INO 

*T.O/ 

1  1 1 
1.11 

Ovenbird 

Seiurus  aurocapillus 

vJ  VCIN 

Migrant 

Mo 
INO 

S  09 

J.7Z 

X 

0  RS 

U.OJ 

Northern  Waterthnish 

Sieurus  noveboracensis 

MOAV  A 
iNvJW  A 

Migrant 

Mo 
INO 

U.JJ 

db 

n  99 
u.zz 

Mourning  Warbler 

Oporomis  Philadelphia 

MOW  A 

Migrant 

Mo 
INO 

n  67 

u.o  / 

-k 

Common  Yellowthroat 

Geothlypis  trichas 

K^yJ  I C 

Migrant 

Mo 
INO 

n  17 

u.l  / 

Z 

nil 

U.  1  i 

Canada  Warbler 

Wilsonia  canadensis 

r^  A\X7  A 

Migrant 

iNO 

U.Uo 

O  OR 

U.UO 

Western  Tanager 

Piranga  ludoviciana 

WETA 

Yes 

2.17 

± 

0.55 

Rose-breasted  Grosbeak 

Pheucticus  ludovicianus 

RBGR 

Migrant 

No 

3.83 

± 

1.04 

Chipping  Sparrow 

Spizella  passerina 

CHSP 

Migrant 

Yes 

5,50 

± 

0.96 

Song  Sparrow 

Melospiza  melodia 

SOSP 

Migrant 

No 

0.17 

± 

0.11 

Lincoln's  Sparrow 

Melospiza  lincolnii 

LISP 

Migrant 

No 

1.83 

± 

0.53 

White-throated  Sparrow 

Zonotrichia  albicollis 

WTSP 

Migrant 

No 

15.92 

± 

1.40 

Dark-eyed  Junco 

Junco  hyemalis 

DEJU 

Migrant 

No 

4.58 

± 

0.81 

Red- winged  Blackbird 

Agelaius  phoeniceus 

RWBL 

Migrant 

No 

0.33 

± 

0.19 

Brown-headed  Cowbird 

Molothrus  ater 

BHCO 

Migrant 

No 

3.08 

± 

0.93 

Pine  Siskin 

Carduelis  pinus 

PISI 

Migrant 

Yes 

1.08 

± 

0.67 
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when  detected  flying  within  the  transect.  This 
species,  unlike  the  other  forest  birds,  does  not 
land  in  trees  or  on  the  ground  while  foraging 
(Godfrey  1986,  Semenchuk  1992).  Other 
species  that  flew  over  during  the  surveys  were 
not  included  in  the  analyses. 

For  each  bird  species,  an  index  of  density  in 
each  1  ha  (2.47  acres)  site  was  calculated  as  the 
maximum  number  of  individuals  detected 
during  the  four  surveys.  Differences  in  bird 
species  richness  and  total  abundance  were 
evaluated  among  stand  groups  using  ANOVA 
for  all  bird  species  together  and  subsequently 
and  for  each  of  the  four  categories  of  bird 
species  (residents  versus  migrants,  and  species 
that  use  big  trees  versus  those  that  do  not). 
Since  bird  species  richness  and  total  abundance 
differed  among  stand  groups  (see  Results), 
differences  in  both  were  evaluated  among 
treatments  after  first  adjusting  for  differences 
among  stand  groups  (Type  III  sum  of  squares  in 
a  2-way  ANOVA;  SAS  Institute  Inc  1989).  In 
all  analyses,  results  were  considered  statistically 
significant  if  the  probability  of  them  occurring 
by  chance  was  less  than  0.05. 

Differences  in  bird  species  richness  and 
abundance  do  not  describe  all  the  differences 
among  bird  communities  because  species 
composition  may  differ  even  when  richness  and 
abundance  are  similar.  I  used  correspondence 
analyses  and  canonical  correspondence  analyses 
(CA  and  CCA  respectively,  Halvorsen  1996)  to 
evaluate  differences  in  bird  communities  among 
stand  groups  and  treatments.  In  all  CAs  and 
CCAs,  the  option  to  down-weight  rare  species 
was  used  (ter  Braak  1992)  because  chance 
events  may  have  resulted  in  these  rare  species 
being  observed  in  non-typical  habitats. 
Differences  among  stand  groups  and  treatments 
were  tested  by  coding  dummy  variables  for  the 
appropriate  categories,  including  these  dummy 
variables  as  environmental  variables  in  the 
analyses,  and  then  testing  for  explained 
variation  using  a  bootstrap  Monte  Carlo  test 
with  99  iterations  (ter  Braak  1992).  To  evaluate 
whether  variation  in  bird  communities  was 
related  to  variation  in  vegetation,  I  included  16 
characteristics  of  live  and  dead  vegetation  in  the 
CCA.  Live  trees  were  categorized  into  three 
groups:  the  dominate  species  (pine  in 


Whitecourt  study  area,  and  aspen  in  Grande 
Prairie  study  area),  the  second  most  dominant 
species  (black  spruce  in  Whitecourt  study  area, 
and  white  spruce  in  Grande  Prairie  smdy  area), 
and  all  other  tree  species.  Within  each  of  these 
three  categories  of  trees,  the  densities  of  trees 
with  DBH  <  and  >  20  cm  (7.9  in.)  was 
calculated.  In  addition,  in  both  study  areas  I 
calculated  densities  of  snags  with  DBH  <  and  > 
20  cm  (7.9  in.),  volume  of  downed  woody 
material  per  ha,  and  percentage  of  the  ground 
surface  covered  by  shrubs  <0.5  m  (<1.6  ft.) 
high,  shrubs  >0.5  m  (>1.6  ft.)  high,  herbs,  grass, 
moss,  and  dead  vegetation. 

Results 

Whitecourt  Study  Area 

Forty-two  species  and  901  individuals  were 
detected  during  the  bird  surveys  (Table  7.1). 
This  bird  community  was  dominated  by  species 
that  migrate  with  only  24%  of  the  species  being 
year-round  residents.  Most  residents  (70%  of  10 
species),  but  fewer  migrants  (40%  of  32 
species)  nested  or  foraged  in  large  live  or  dead 
trees.  More  bird  species  (Table  7.2)  and 
individuals  (Table  7.3)  were  observed  in  stand 
group  B  than  in  the  other  three  stand  groups. 
That  difference  was  largely  due  to  more  resident 
bird  species  that  use  large  trees  being  present  in 
stand  group  B.  Species  richness  and  abundance 
were  similar  among  treatments  for  all  categories 
of  birds  after  adjusting  for  differences  among 
stand  groups  (Tables  7.2  and  7.3). 

When  evaluating  differences  in  bird 
communities  sunong  sites,  the  first  and  second 
functions  of  the  CA  had  eigen  values  of  0.07 
and  0.04,  respectively.  Combined,  these  two 
functions  accounted  for  only  3.4%  of  the 
variation  in  bird  species  density.  However,  even 
with  this  small  amount  of  explained  variation, 
bird  communities  were  statistically  different 
among  stand  groups  (Fig.  7.1;  F  =  1.8,  P  < 
0.001).  After  adjusting  for  differences  among 
stand  groups,  bird  communities  did  not  differ 
among  treatments  (during  the  bootstrap  Monte 
Carlo  test  permutations  were  restricted  to  within 
the  appropriate  stand  group;  F  =  1.2,  P  =  0.10). 
In  addition,  all  categories  of  bird  species  were 
present  in  all  stand  groups  (Fig.  7.2).  Variation 
in  the  bird  communities  was  not  related  to 
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TABLE  7.2.  Mean  (±  SE)  species  richness  within  sites  within  each  stand  group  and  treatment. 
Species  were  categorized  by  whether  or  not  they  used  large  trees  or  snags  and  their  migration 
habitats.  There  were  3  and  121  degrees  of  freedom  for  the  tests  involving  stand  groups  and  2  and 
1 19  degrees  of  freedom  for  the  tests  involving  treatments.  Groups  with  different  superscripts  are 
statistically  different  (Student-Neuman-Keuls  range  test). 


Categories  Use  Big  Trees  Do  Not  Use  Big  Trees 

Test   Total 

Statistic      Resident        Migrant         Resident  Migrant 


Whitecourt  Study  Area 

Stand  Group  A  0.2  ±0.1'' 

Stand  Group  B  0.5  ±0.1* 

Stand  Group  C  0.2  ±0.1** 

Stand  Group  D  0.2  ±0.1'' 

F  3.5 

P  0.02 

Treatment  1  0.2  ±0.1 

Treatment  2  0.3  ±0.1 

Treatments  0.3  ±0.1 

F  0.6 

P  0.94 

Grande  Prairie  Study  Area 

Stand  Group  A  0.1  ±0.1 

Stand  Group  B  0.1  ±0.1 

Stand  Group  C  0.2  ±0.1 

Stand  Group  D  0.1  ±0.1 

F  0.2 

P  0.92 

Treatment  1  0.1  ±0.1 

Treatment  2  0.1  ±0.1 

Treatments  0.1  ±0.1 

F  0.6 

P  0.94 


0  Q    n  0 

Z.7  3  u.z 

n  A  -4-  n  1 
u.o  X  u.  1 

1  1  +  n  1 

i .  1  X  U.  1 

*t.o  X  U.J 

-3  rj  +  n  9 

u.o  X  w.  1 

i.j  X  u.z 

J.O  X  U.*T 

2.8  ±0.2 

0.7  ±0.1 

1.5  ±0.2 

5.1  ±0.3"^ 

2.3  ±0.3 

0.6  ±0.1 

1.4  ±0.2 

4.5  ±0.4** 

l.O 

1  ft 

l.O 

n  Q 

9  AQ 

u.zz 

0  AQ 

U.T-l/ 

U.UJ 

Z.J  Jt  u.z 

kJ.o  z:  U.  1 

1  .z  X  u.z 

M-.  /  X  U.J 

0  o    n  9 

KJ.D  3:  U.  1 

1  *;  -»-  n  9 
i.j  X  u.z 

J.J  X  U.J 

2.8  ±0.2 

0.6  ±0.1 

1.4  ±0.2 

5.0  ±0.3 

0.8 

0.3 

1.2 

1.3 

0.43 

0.77 

0.32 

0.27 

3.2  ±0.2 

0.3  ±0.1'** 

4.1  ±0.2** 

7.7  ±0.4" 

2.7  ±  0.2 

0.4  ±0.1' 

4.0  ±0.2'' 

7.3  ±0.4" 

3.4  ±0.2 

0.1  ±0.1" 

5.3  ±0.3" 

9.0  ±0.5* 

2.8  ±  0.2 

0.3  ±0.1'* 

5.0  ±0.3" 

8.2  ±0.4* 

2.5 

3.2 

6.4 

4.6 

0.06 

0.03 

0.001 

0.001 

3.0  ±0.2 

0.3  ±0.1 

4.5  ±  0.2 

7.9  ±0.04 

3.2  ±0.2 

0.2  ±0.1 

5.0  ±0.3 

8.5  ±0.06 

2.9  ±  0.2 

0.3  ±0.1 

4.3  ±  0.3 

7.6  ±0.04 

0.8 

0.3 

2.3 

2.2 

0.44 

0.77 

0.10 

0.14 
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TABLE  7.3.  Mean  (±  SE)  abundance  per  ha  within  each  stand  group  and  treatment.  Species 
were  categorized  by  whether  or  not  they  used  large  trees  or  snags  and  their  migration  habitats. 
There  were  3  and  121  degrees  of  freedom  for  the  tests  involving  stand  groups  and  2  and  1 19 
degrees  of  freedom  for  the  tests  involving  treatments.  Groups  with  different  superscripts  are 
statistically  different  (Student-Neuman-Keuls  range  test). 


Categories 

Use  Big  Trees 

Do  Not  Use  Big  Trees 

Test 

Total 

Statistic 

Resident 

Migrant 

Resident 

Migrant 

Whitecourt  Study  Area 

Stand  Group  A 

0.4  ±  0.2 

3.3  ±0.2 

0.6  ±0.1 

1.2  ±0.1 

5.5  ±0.4 

Stand  Group  B 

0.6  ±0.1 

3.5  ±  0.3 

0.7  ±0.1 

1.6  ±0.2 

6.3  ±0.5 

Stand  Group  C 

0.3  ±  0.2 

3.2  ±0.3 

0.7  ±0.1 

1.6  ±0.2 

5.9  ±0.4 

Stand  Group  D 

0.2  ±0.1 

2.7  ±  0.3 

0.7  ±0.1 

1.5  ±0.2 

5.1  ±0.4 

F 

1.2 

1.8 

0.1 

0.8 

2.3 

P 

0.31 

0.14 

0.95 

0.46 

0.08 

Treatment  1 

0.3  ±0.1 

2.9  ±  0.2 

0.8  ±0.1 

1.2±0.1 

5.2  ±0.4 

Treatment  2 

0.5  ±0.1 

3.4  ±  0.3 

0.5  ±0.1 

1.7  ±0.2 

6.2  ±0.5 

Treatment  3 

0.3  ±0.1 

3.2  ±  0.2 

0.7  ±0.1 

1.4  ±0.2 

5.6  ±0.4 

F 

0.4 

1.3 

0.5 

1.7 

2.0 

P 

0.65 

0.29 

0.61 

0.19 

0.11 

Grande  Prairie  Study  Area 

Stand  Group  A 

0.1  ±0.1 

3.9  ±0.3*" 

0.3  ±0.1*" 

5.0  ±0.3' 

9.3  ±0.5" 

Stand  Group  B 

0.1  ±0.1 

3.4  ±  0.3 

0.6  ±0.2* 

4.7  ±0.3' 

8.8  ±0.4" 

Stand  Group  C 

0.2  ±0.1 

4.4  ±0.4* 

0.1  ±0.1" 

7.4  ±0.5* 

12.0  ±0.6* 

Stand  Group  D 

0.2  ±0.1 

3.1  ±0.2" 

0.3  ±0.1*" 

6.2  ±0.4" 

9.8  ±0.5" 

F 

0.3 

3.3 

3.2 

11.8 

9.7 

P 

0.86 

0.02 

0.03 

<0.001 

<0.001 

Treatment  1 

0.1  ±0.1 

3.5  ±  0.2 

0.3  ±0.1 

5.6  ±0.2 

9.6  ±0.5 

Treatment  2 

0.1  ±0.1 

4.3  ±0.4 

0.3  ±0.1 

6.2  ±0.3 

10.7  ±0.6 

Treatment  3 

0.2  ±0.1 

3.3  ±0.3 

0.3  ±0.1 

5.5  ±0.4 

9.6  ±0.5 

F 

0.4 

25 

0.1 

1.5 

2.0 

P 

0.66 

0.09 

0.87 

0.22 

0.11 
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FIGURE  7.1.  Correspondence  analyses  based  on  bird  densities  in  the  Whitecourt  study  area 
during  1997.  Sites  that  are  close  together  within  the  correspondence  plot  have  similar  bird 
communities  and  sites  that  are  far  apart  have  dissimilar  bird  communities.  Sites  were  categorized 
based  on  stand  group.  For  each  stand  group  bivariate  ellipses  that  enclose  90%  of  the  sites  were 
presented  (Systat  Inc.  1989). 
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3i 


^    B     Use  Big  Trees 
on     Do  Not  Use  Big  Trees 
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□ 
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□ 
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□ 
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□ 
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□ 
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□ 
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□ 


-1  0  1 

Correspondence  Function  1 


FIGURE  7.2.  Predicted  location  for  each  bird  species  based  on  the  correspondence  analyses  of 
bird  densities  in  the  Whitecourt  study  area  during  1997.  Species  that  are  close  together  have 
strong  positive  covariance  in  abundance  among  sites  whereas  species  that  are  far  apart  have 
strong  negative  covariance  in  abundance  among  sites.  Species  were  categorized  based  on  their 
migration  and  nesting/foraging  habits. 
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variation  in  vegetation  characteristics  among 
sites  (F=  1.1,  P  =  0.32). 

Grande  Prairie  Study  Area 

Fifty  species  and  1592  individuals  were  detected 
during  the  bird  surveys  (Table  7.1).  This  bird 
community  was  dominated  by  species  that 
migrate  with  only  18%  of  the  species  being 
year-round  residents.  Most  residents  (67%  of  9 
species),  but  fewer  migrants  (34%  of  41 
species)  nested  or  foraged  in  large  live  or  dead 
trees.  Both  species  richness  (Table  7.2)  and 
number  of  individuals  (Table  7.3)  differed 
among  stand  groups.  Differences  in  species 
richness  were  mainly  due  to  differences  for 
birds  that  did  not  use  large  trees  and  snags, 
whereas  differences  in  abundance  were  due  to 
differences  for  birds  that  did  not  use  large  trees 
and  snags  and  differences  for  birds  that  migrate 
and  use  large  trees.  After  adjusting  for 
differences  in  richness  and  abundance  among 
stand  groups,  richness  and  abundance  were 
similcu-  among  treatments  for  all  categories  of 
birds. 

When  evaluating  differences  in  bird 
communities  among  sites,  the  first  and  second 
functions  of  the  CA  had  eigen  values  of  0.07 
and  0.04,  respectively.  Combined,  these  two 
functions  accounted  for  only  4.1%  of  the 
variation  in  bird  species  density.  However,  even 
with  this  small  amount  of  explained  variation, 
bird  communities  were  statistically  different 
among  stand  groups  (Fig.  7.3;  F  =  1 .9,  P  < 
0.01).  After  adjusting  for  differences  among 
stand  groups  bird  communities  did  not  differ 
among  treatments  (F  =  1.2,  P  =  0.14).  In 
addition,  all  categories  of  bird  species  were 
present  in  all  stand  groups  (Fig.  7.4).  Variation 
in  the  bird  communities  was  not  related  to 
variation  in  vegetation  characteristics  among 
sites  (F=  1.2,  P  =  0.18). 

Discussion 

To  achieve  a  rigorous  experimental  design,  we 
chose  two  study  areas  that  were  large  and 
relatively  homogenous.  Within  these  study 
areas,  the  stands  (i.e.,  experimental  harvest 
units)  were  similar  to  each  other  prior  to 
harvest,  but  as  expected  within  natural  forests, 
the  stands  were  not  identical  to  each  other.  We 


created  four  groups  of  stands  in  each  study  area 
based  on  vegetation  in  and  around  the  stands 
(see  Chapter  2).  This  resulted  in  a  study  design 
where  some  of  the  variation  in  vegetation  and 
structure  could  be  accounted  for  statistically  in  a 
stratified  analysis  (Chapters  3  and  4).  Similarly, 
within  both  study  areas,  variation  among  bird 
communities  was  related  to  differences  among 
stand  groups.  That  was  true  for  both  the  simple 
analyses  of  bird  species  richness  and  abundance, 
and  for  the  more  complex  community  analyses. 
However,  variation  among  stand  groups 
accounted  for  only  a  small  amount  of  the 
variation  in  bird  communities  and  the  bird 
communities  overlapped  greatly  among  the 
stand  groups.  This  was  expected  since  there 
were  only  slight  differences  in  vegetation 
among  stand  groups. 

After  controlling  for  the  differences  among 
stand  groups,  bird  communities  did  not  differ 
among  treatments  prior  to  harvest. 
Consequently,  if  bird  communities  differ  among 
treatments  after  harvest,  those  differences  can 
be  attributed  unambiguously  to  differences  in 
the  size  and  dispersion  of  tree  patches  that  were 
retained  in  these  experimental  areas  during 
harvest.  Thus,  as  planned,  this  study  will 
provide  a  rigorous  comparison  of  bird 
communities  in  among  the  three  dispersions  of 
standing  trees  and  snags.  The  bird  communities 
found  in  the  stands  prior  to  harvest  are 
benchmarks  from  which  I  will  evaluate 
differences  after  harvest. 

Future  Research 

Bird  communities  will  be  surveyed  during  1998 
and  1999.  However,  most  of  the  post-harvest 
bird  surveys  will  be  conducted  during  1999 
because  birds  may  crowd  into  the  remaining 
unharvested  habitat  immediately  after  harvest 
and  those  high  densities  probably  will  decrease 
during  the  next  year  (Lovejoy  et  al.  1986, 
Schmiegelow  1997).  Bird  surveys  during  1999 
will  be  an  extended  version  of  those  conducted 
prior  to  harvest.  Within  each  stand,  the  10  ha 
(24.7  acres)  study  area  (Figs.  2.1  and  2.2)  will 
be  surveyed  using  a  combination  of  strip 
transects  and  point  counts  (see  Methods).  In 
addition,  eight  point  counts  will  be  centered  in 
residual 
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FIGURE  7.3.  Correspondence  analyses  based  on  bird  densities  in  the  Grand  Prairie  study  area 
during  1997.  Sites  that  are  close  together  within  the  correspondence  plot  have  similar  bird 
communities  and  sites  that  are  far  apart  have  dis-similar  bird  communities.  Sites  were 
categorized  based  on  stand  group.  For  each  stand  group  bivariate  ellipses  that  enclosed  90%  of 
the  sites  were  presented. 
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FIGURE  7.4.  Predicted  location  for  each  bird  species  based  on  the  correspondence  analyses  of 
bird  densities  in  the  Grande  Prairie  study  area  during  1997.  Species  that  are  close  together  have 
strong  positive  covariance  in  abundance  among  sites  whereas  species  that  are  far  apart  have 
strong  negative  covariance  in  abundance  among  sites.  Species  were  categorized  based  on  their 
migration  and  nesting/foraging  habits. 
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clumps  within  each  stand  and  four  counts  will 
be  located  >100  m  (>109  yd.)  from  residual 
clumps  in  each  stand.  .Surveys  will  be  conducted 
in  all  treatments  once  between  15  and  30  April, 
and  three  times  at  bi-weekly  intervals  between 
15  May  and  30  June.  During  the  reduced  1998 
survey,  only  40  point  counts  will  be  conducted 
in  each  study  area.  Ten  of  these  point  counts 
will  be  centered  within  small  residual  patches, 
ten  centered  within  mid-sized  residual  patches, 
ten  centered  within  large  residual  patches,  and 
ten  >100  m  (>109  yd.)  from  residual  patches. 
Results  from  these  bird  surveys  will  be  analyzed 
and  presented  during  fall  and  winters  of  1998 
and  1999. 
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CHAPTER  8.  THE  USE  OF  RESIDUAL  PATCHES  BY  BATS  IN  THE  BOREAL 
FOREST  OF  NORTHWESTERN  ALBERTA 


Lui  Marinelli 


introduction 

Timber  harvesting  strategies  are  restructuring 
the  boreal  forest  mosaic  in  Alberta.  With  the 
recent  introduction  of  aspen  into  harvesting 
plans,  forests  are  being  reshaped  at  an  even 
faster  rate.  Rolling  (1995)  characterized 
traditional  resource  management  as  controlling 
the  variability  of  ecosystems  so  that  a  steady 
flow  of  goods  and  services  could  be  produced 
for  humans.  This  traditional  view  of  forest  value 
is  limited  to  the  monetary  worth  of  present  and 
future  timber  resources.  Emerging 
environmental,  ecological,  and  social  concerns 
are  forcing  a  much  broader  perspective  of 
resource-use  strategies.  Strategies  of  ecosystem 
management  are  being  developed  (e.g..  Alberta 
Forest  Conservation  Strategy  1997), 
encouraging  private  companies  to  reassess  the 
value  of  the  forest  and  modify  traditional 
harvesting  practices  (e.g..  Alberta  Pacific  Forest 
Industries  Inc.  1996). 

Harvesting  practices  are  evolving  to  include 
leaving  standing  trees  and  snags  in  harvested 
areas  in  an  effort  to  retain  pre-harvest  species  of 
flora  and  fauna.  Such  forest  structure  may  make 
harvested  areas  more  desirable  to  native  bats 
that  use  the  forest  environment  to  forage,  roost, 
and  breed  (Thomas  1988).  Currently,  the 
pattern  in  which  residual  trees  and  snags  are  left 
on  the  landscape  is  decided"  subjectively 
(Alberta  Pacific  Forest  Industries  Inc.  1996). 
Objective  decision-making  requires  a  greater 
understanding  of  the  effect  on  biotic 
communities  of  retaining  residual  elements  and 
the  manner  of  their  dispersion. 

The  boreal  forest  of  Alberta  supports  five 
species  of  bats  (van  Zyll  de  Jong  1983).  These 
bats  show  a  fair  degree  of  fidelity  to  roost  sites 
in  the  forest  when  they  return  from 
hibemaculum  or  migration.  During  the  spring 
and  summer,  cavity-roosting  colonial  species 
use  a  wide  variety  of  cavities  including  spaces 


behind  exfoliated  bark  and  natural  or  excavated 
cavities  in  large,  old  snags  (Thomas  and  La  Val 
1988).  Within  Douglas-fir  forests  of  southern 
Washington  State,  old  forests  were  considered 
critical  for  bats  because  they  contained  large, 
decayed  snags,  which  provide  crevices  for 
cavity-roosting  bat  species  (Thomas  1988).  The 
two  foliage-roosting  solitary  species  {Lasiurus 
cinereus  and  Lasionycteris  noctivagans) 
generally  do  not  exploit  cavities  but  roost  in 
fixed  or  shifting  locations  in  foliage  of  large, 
live  trees  (Barclay  1984).  In  some  species, 
males  and  females  use  different  habitats 
(Barclay  1991).  Males  do  not  maintain  high 
body  temperatures  during  the  day  and  may 
select  cool  roosts  to  minimize  their  metabolic 
rates  and  diurnal  energy  expenditures  (Thomas 
1988).  Females  have  a  reproductive  constraint 
on  their  ability  to  exploit  the  energy  saving  of 
torpor  because  fetal  growth  rates  and 
presumably  milk  production  are  temperature- 
dependent  processes  (Racey  1982).  During 
pregnancy  and  lactation,  females  maintain  body 
temperatures  well  above  ambient  during  the  day 
to  maximize  fetal  growth  rates  (Kunz  1982b, 
Racey  1982)  and  as  such  tend  to  select  warm 
roost  sites.  Thus,  retaining  bats  in  forest  patches 
will  depend  upon  requirements  of  different  bat 
species  and  sexes. 

In  this  study,  three  treatments  with  different 
clump  size  and  distances  between  clumps  in  two 
experimentally  harvested  areas  will  be 
established  to  evaluate  how  bat  communities 
respond  to  harvested  forests.  We  will 
concentrate  on  differences  that  are  detected  in 
the  first  two  years  following  harvest,  however, 
throughout  the  forest  rotation  the  opportunity 
exists  to  evaluate  the  value  to  bats  of  the  patches 
that  are  created  in  the  experiment. 

Objectives 

I  will  compare  bat  communities  among  three 
different  harvesting  treatments  (Chapter  2)  to 
evaluate  which  pattern  of  dispersion  of  residual 
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trees  maximizes  bat  diversity  and  abundance. 
Pre-harvest  sampling  will  be  used  to  assess 
whether  bat  communities  were  similar  among 
treatments  prior  to  harvest.  Post-harvesting  bat 
communities  will  be  compared  among  the  three 
treatments  using  data  that  will  be  collected 
during  1998  and  1999. 

Methods 

A  description  of  the  study  areas  and 
experimental  harvesting  is  presented  in  the 
Study  Area  and  Methods  section  (Chapter  2). 
Descriptions  of  the  pre-harvest  trees,  snags, 
downed  woody  materials,  and  understory 
vegetation  cu-e  presented  in  Chapters  3-5. 

Bat  detectors  were  used  to  estimate  species 
richness  and  abundance  in  each  stand  in  both 
Whitecourt  and  Grande  Prairie  smdy  areas  in 
July  and  August  1997.  Ultrasonic  detectors  do 
not  provide  data  that  can  be  translated  directly 
into  estimates  of  population  density  (Thomas 
and  La  Val  1988);  they  do  however,  provide  a 
relatively  unbiased  index  of  the  levels  of  use  or 
animal  abundance  among  habitats  (Thomas 
1988).  Detectors  consisted  of  an  Anabat  IV 
ultrasonic  bat  detector  (Titley  Electronics, 
AUS),  a  delay  switch,  a  cassette  recorder,  and  a 
12- volt  battery,  all  of  which  were  housed  in  a 
1.5  litre  (1.6  gal.)  rubbermaid  container.  A  2  cm 
(0.8  in.)  hole  was  cut  into  the  side  of  the 
container  to  accommodate  the  microphone.  The 
battery  was  used  to  run  the  detector  and  the 
delay  switch.  The  delay  switch  was  set  to  turn 
the  detector  on  at  dusk  and  off  at  dawn.  The 
sensitivity  setting  on  the  detectors  was  set 
between  6  and  7  to  minimize  insect  noise  and  to 
eliminate  detection  of  bats  flying  in  adjacent 
stand  groups.  To  allow  the  bat  call  to  be 
recorded  by  a  cassette  recorder,  the  frequency  of 
each  call  was  reduced  from  the  ultrasonic  to  the 
audible  range  by  setting  the  division  ratio  on 
detectors  to  16. 

Bats  were  surveyed  from  three  points, 
approximately  300  m  (328  yd.)  apart,  in  each 
stand.  Two  stands  in  each  study  area  were 
sampled  simulteineously.  Detectors  were 
attached  to  trees  1.5  m  (5  ft.)  off  the  ground  and 
oriented  such  that  the  microphone  was  pointing 
45-degree  above  horizontal.  Detectors  were  run 


for  three  consecutive  nights  in  each  stand  and 
were  checked  each  morning  to  ensure  proper 
functioning.  After  three  days,  detectors  were 
moved  to  new  stands. 

Bats  were  identified  by  analyzing  their  search- 
phase  echolocation  calls.  Each  sequence  of  one 
or  more  echolocation  pulses  with  <1  s  between 
sequential  pulses  was  recorded  as  a  pass  by  a 
bat  (Fenton  1970).  Calls  were  analyzed  using 
Anabat  cind  Analook  softwcire  (Titley 
Electronics,  AUS).  Frequency  range  and  call 
shape  from  the  calls  we  recorded  were 
compared  with  a  library  of  calls  from  known 
species. 

Differences  in  bat  abundance  were  evaluated 
among  stand  groups  and  treatments  using  one- 
way ANOVAs.  In  all  analyses,  results  were 
considered  statistically  significant  if  the 
probability  of  them  occurring  by  chance  was 
less  than  0.05.  Canonical  correlation  analysis 
was  used  to  evaluate  the  relationship  between 
the  bat  abundance  and  characteristics  of  trees  in 
each  stand.  In  both  areas,  I  included  the 
following  tree  variables  in  the  analyses:  density 
of  snags  with  a  diameter  at  breast  height  (DBH) 
of  <20  cm  (<7.9  in.),  snags  >20  cm  (>7.9  in.) 
DBH,  and  total  volume  of  down  woody 
material.  In  the  Whitecourt  study  area,  density 
of  pine  and  black  spruce  <20  cm  (<7.9  in.)  DBH 
and  density  of  pine  and  black  spruce  >20  cm 
(>7.9  in.)  DBH  were  included  in  the  analyses. 
In  the  Grande  Prairie  study  area,  density  of 
aspen  and  white  spruce  <20  cm  (<7.9  in.)  DBH 
and  density  of  aspen  and  white  spruce  >20  cm 
(>7.9  in.)  DBH  were  included  in  the  analyses. 

Results 

Based  on  756  hrs  of  data  collection  in  each  area 
almost  twice  as  many  bat  calls  were  detected  in 
the  Grande  Prairie  study  area  (n  =  120)  as  in  the 
Whitecourt  study  area  (n  =  67).  Of  the  bats 
detected  (n  =  187),  the  majority  were  identified 
as  little  brown  bats  {Myotis  lucifugus).  Hoary 
bats  (Lasiurus  cinereus)  were  detected  three 
times  in  the  Whitecourt  study  area. 
Accordingly,  I  restricted  the  analyses  to 
detections  of  the  little  brown  bats. 


03/10/98 


Client  copy  -  not  for  citation  or  reproduction 


68 


Use  of  Residual  Patches  by  Bats 


TABLE  8.1.  Mean  (±  SE)  number  of  little  brown  bat  calls  detected  within  each  stand  group  and 
treatment  in  Whitecourt  and  Grande  Prairie  study  areas. 


Whitecourt 

Grande  Prairie 

Stand  group  1 

2.3  (1.5) 

14.0  (4.4) 

2 

2.7  (1.5) 

7.3  (2.0) 

3 

2.7(1.2) 

1.7  (0.7) 

4 

8.0  (6.5) 

2.3  (1.2) 

Treatment  1 

1.8(1.1) 

6.0  (2.3) 

2 

2.3(1.0) 

7.5  (5.1) 

3 

7.8  (4.5) 

5.5  (2.2) 
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Bat  detections  in  the  Whitecourt  study  area  did 
not  differ  among  treatments  (F^,,  =  1.50,  P  = 
0.27)  or  among  stand  groups  (F, =  0.62,  P  = 
0.62;  Table  8.1).  In  the  Grande  Prairie  study 
area  no  significant  difference  was  detected 
among  treatments  (F^  „  =  0.09,  P  =  0.92). 
However,  among  stand  groups,  some 
differences  were  significant  (F3  „  =  5.23,  P  = 
0.03).  A  Tukey  post-hoc  test  suggested  that 
difference  between  stand  group  1  and  stand 
groups  3  and  4  are  responsible  for  the 
significant  difference. 

The  density  of  trees  and  snags  accounted  for 
only  a  small  amount  of  the  variation  in  the  bat 
numbers  in  both  study  areas  (4%  in  Grande 
Prairie,  and  5%  in  Whitecourt).  In  the  Grande 
Prairie  study  area,  aspen  density  [both  <  and  > 
20  cm  (7.9  in.)  DBH]  was  the  best  predictor  of 
bat  abundance  (Fig.  8.1).  Snag  density  [>20  cm 
(>7.9  in.)  DBH]  v/as  the  best  predictor  of  bat 
abundance  in  the  Whitecourt  study  area  and  was 
positively  correlated  with  bat  abundance  (Fig. 
8.2).  The  graphical  depiction  of  the  relationship 
between  bat  abundcince  and  tree  veiriables  (Fig. 
8.1  and  8.2)  also  shows  the  lack  of 
differentiation  among  the  stand  groups  and 
treatments  in  both  study  areas. 

Discussion 

Stands  in  each  area  were  categorized  into  four 
stand  groups  based  on  observed  similarity  in 
vegetation  in  and  around  the  stands  (see  Chapter 
2).  Little  brown  bats  were  equally  abundant  in 
all  stand  groups  in  the  Whitecourt  area,  but  were 
more  abundant  in  stand  group  1  in  the  Grande 
Prairie  area.  Accordingly,  if  the  density  of  little 
brown  bats  in  the  Whitecourt  area  differs  among 
stands  after  harvest,  those  differences  can  be 
attributed  to  differences  in  the  size  and 
dispersion  of  tree  patches  that  were  retained  at 
harvest.  However,  in  the  Grande  Prairie  study 
area,  post-harvest  differences  will  have  to  be 
evaluated  after  taking  into  account  pre-harvest 
differences. 

Aspen  trees  in  the  Grande  Prairie  study  area  and 
snags  in  the  Whitecourt  study  area  appeared  to 
be  important  in  predicting  bat  abundance.  Bats 
rely  on  crevices  in  trees  or  raised  bark  to  roost 
during  the  day.  Accordingly,  snags  in  the 


Whitecourt  study  area  may  be  important 
because  of  the  roost  sites  they  offer.  Bats  can 
also  roost  in  live  trees,  using  naturally  occurring 
cracks  or  excavated  cavities.  Bats  may  also 
select  an  area  based  on  the  dispersion  of  live 
trees.  Although  bats  do  not  appear  to  spend  a 
significant  portion  of  their  time  feeding  in  the 
forest  (Thomas  1988),  when  they  do,  they 
require  a  corridor  to  fly  and  feed.  Mature  aspen 
in  the  Grande  Prairie  study  area  may  provide 
roost  sites  and/or  the  proper  corridors  for 
feeding. 

Future  Research 

Bat  communities  will  be  surveyed  during  1998 
and  1999  to  evaluate  the  effect  of  the  harvest 
treatments  on  bats.  Three  bat  detectors  units 
will  be  positioned  in  a  stand  each  night;  one  at 
the  center  of  a  residual  patch,  one  at  the  patch 
edge,  and  one  within  the  clear-cut  area.  Bat 
calls  will  be  recorded  for  three  consecutive 
nights.  Flying  insect  communities  will  be 
sampled  concurrently  at  these  same  locations 
using  insect  intercept  traps  that  will  be 
suspended  at  a  height  of  2  m  (6.6  ft.).  Insect 
collection  will  be  limited  to  the  time  when  bats 
are  most  active,  from  sunset  to  three  hours  after 
sunset,  by  using  automated  traps  that  will  be 
controlled  by  photocells,  timers,  and  motors. 
Results  from  bat  surveys  will  be  analyzed  and 
presented  during  fall  and  winters  of  1998  and 
1999. 
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FIGURE  8.1.  Canonical  correlation  analyses  based  on  abundance  of  little  brown  bats  in  boreal 
forest  stands  near  Grande  Prairie,  AB  during  1997.  The  first  correlation  function  describes  an 
increasing  density  of  aspen  trees.  The  symbols  represent  each  stand  identified  by  stand  group 
(A-D)  and  treatment  (1-3). 


03/10/98 


Client  copy  -  not  for  citation  or  reproduction 


71 


Use  of  Residual  Patches  by  Bats 


LU 
O 

z 
< 
o 
z 

m 
< 

< 
ffl 

LU 
> 


LU 


2.00  - 


1.00 


0.00 


•1.00 


D3 

83 

^«  ni 

C3     ^ ' 
A2D1 

B2 
B1 

D2  A3 
C2 

A1 

•1.50 


■0.50  0.50  1.50 

CORRELATION  FUNCTION  1 


2.50 


FIGURE  8.2.  Canonical  correlation  analyses  based  on  abundance  of  little  brown  bats  in  boreal 
forest  stands  near  Whitecourt,  AB  during  1997.  The  first  correlation  function  describes  an 
increasing  density  of  snags  with  a  DBH  of  >20  cm  (>7.9  in.).  The  symbols  represent  each  stand 
identified  by  stand  group  (A-D)  and  treatment  (1-3). 
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